arXiv: 1501.04095v2 [astro-ph.CO] 3 Apr 2015 


Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 7 April 2015 (MN IAT^K style file v2.2) 


The galaxy cluster outskirts probed by Chandra 


Andrea Morandi 1 *, Ming Sun 1 , William Forman 2 , Christine Jones 2 

1 Physics Department, University of Alabama in Huntsville, Huntsville, AL 35899, USA 

2 Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138, USA 


7 April 2015 


ABSTRACT 

We studied the physical properties of the intracluster medium in the virialization re¬ 
gion of a sample of 320 clusters (0.056 < 2 < 1.24, kT )> 3 keV) in the Chandra 
archive. With the emission measure profiles from this large sample, the typical gas 
density, gas slope and gas fraction can be constrained out to and beyond R 200 ■ We ob¬ 
serve a steepening of the density profiles beyond i?soo with /3 ~ 0.68 at i? 50 o and /3 ~ 1 
at i ?200 and beyond. By tracking the direction of the cosmic filaments approximately 
with the ICM eccentricity, we report that galaxy clusters deviate from spherical sym¬ 
metry, with only small differences between relaxed and disturbed systems. We also 
did not find evolution of the gas density with redshift, confirming its self-similar evo¬ 
lution. The value of the baryon fraction reaches the cosmic value at i? 200 : however, 
systematics due to non-thermal pressure support and dumpiness might enhance the 
measured gas fraction, leading to an actual deficit of the baryon budget with respect 
to the primordial value. This study has important implications for understanding the 
ICM physics in the outskirts. 
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1 INTRODUCTION 

Exploring the virialization region of galaxy clusters is vi¬ 
tal to our understanding of large scale structure formation, 
offering a direct view of cluster growth. This topic has re¬ 
cently raised the attention of the scientific community, since 
an accurate measurement of the mass and baryonic fraction 
of cluster outskirts can provide constraints on cosmological 
parameters, an d allow us to test th e validity of the CDM 
framework (e.g. iReiprich et al.ll2013h . 

Nevertheless, the outskirts remain a relatively unex¬ 
plored territory, bearing the signature of a complex physics 
such as ongoing accretion processes, significant d epartures 
from virialization and hy drostatic equilibrium (ILau et al.l 
l2009h . clumping of the gas llMorandi fc Cuill2014h . These ef¬ 
fects may bias, in turn, the measurements of cluster masses, 
thus limiting the use of galaxy clusters as a high-precision 
cosmological proxy. Hence a deeper understanding of the 
state of the intracluster gas in cluster outskirts is required. 

Observations are very challenging in cluster outskirts, 
since the X-ray surface brightness drops below the back¬ 
ground level at large radii. Thanks to its low particle back¬ 
ground from its low orbit, Suzaku observations have ex¬ 
tended X-ray measurements of the intracluster medium 
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(ICM) profiles out to and beyond R 200 , where R 200 is the 
radius within which the mean total density is 200 times the 
critical density of the Universe. Initial results were surpris¬ 
ing, revealing significant departures from the theoretical pre¬ 
dictions. These observations showed that, at large cluster¬ 
centric radii (R )> R 200 ), the observ e d entropy profiles 
becom es flatter (e.g., iBautz et all 120091 iKawaharada et aD 
l2010h . and the observed gas fract ion exceeded the cosmi c 
baryon fraction of the universe (ISimionescu et ahl 1201 il l. 
Suzaku measurements call for appreciable gas dumpiness in 
the outs kirts, infe rred from its effects on thermodynamic 
profiles dWalker et ahlEoU l or from the high values, larger 
than the cosmic baryon fraction, of the gas mass fraction 
towards the virial radius. However, difficulties arise in mod¬ 
eling out point sources and/or galactic foreground with 
the limited angular resolution (~ 2 arcmin) of Suzaku at 
the level required for measuring the extremely low surface 
brightness in the cluster outskirts. In particular, the Suzaku 
point-spread function and the stray light can cause signif¬ 
icant contamination from the brightness in the central re¬ 
gions to large radii. Moreover, Suzaku observations have 
been mostly performed along narrows arms, which might re¬ 
flects preferential directions connected with the large-scale 
structure (e.g., in the direction of filaments). Therefore, it 
is possible that the aforementioned measurements are not 
representative of the cluster as a whole. 
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Independent ROSAT PSPC observations (with low 
background and large field of view) have shown the steep¬ 
ening of the gas density profile at large radii, but at the 
level considerably larger than those inferred from the Suzaku 
observations (the ICM slope (3 = 0.890 ± 0.026 at R200, 
lEckert et al.ll2012h . 


I 11 this respect, independent measurements with 
Chandra with respect to Suzaku can shed light on 
syste matic uncertainties in th e current measur e ments 
Ettori fc Balestral l2009bl : lEttori fc MolendH l201ll : 
iMoretti et al.ll2011f ): indeed, despite its higher particle back¬ 
ground, Chandra provides a superior angular resolution 
which is fundamental in order to: i) remove emission from 
unrelated sources, e.g. point sources; ii) to constrain the 
emission of clusters out to the virial radius, especially for 
higher-redshift cool-core clusters, for which there might be 
a non-negligible contribution from the bright cluster core 
to the emission in the outer volumes and from secondary 
scatter by sources outside the field of view. Note that the 
properties of the particle background are well known via the 
stowed background observations, with uncertainties in the 
backg round across its spectrum w hich are within 2 per¬ 
cent flHickox fc Markevitchl f20061 ') and see also discussion in 
ISun et alJ d2009h~ 


This data analysis project represents a follow-up of our 
previous works on studying cluster thermodynamic proper¬ 
ties, including emiss ion measure, gas dens i ty, and gas frac¬ 
tion, out to R200 dMorandi et al. l2013bl : (Morandi fc Cuil 
120141 '). Indeed, when data of sufficient depth are avail¬ 
able, we have shown that Chandra can accurately mea¬ 
sure the sur face brightness and temperature of clusters 
out to Rr >no dBo nament e et al.l 120131 : IMorandi et aid l2013bl : 
iMorandi fc Cuill2014l . e.g.. in Abell 1835 and Abell 133). 
The Chandra’s superior angular resolution enables robust 
identification and removal of point sources from the X-ray 
images, while minimizing the contribution from the bright 
cluster core to the emission in the outer volumes - both 
challenging tasks for Suzaku - and to provide an inde¬ 
pendent confirmation of the Suzaku results. Observational 
constraints on the gas clumping factor inferred from the 
inhomogeneities of the X-ray surface brightness are also 
significantly smaller than those inferr ed from the Suzaku 
measurements (C ~ 1.5 — 2 at R 200 , Eckert et al.ll2013al : 
IMorandi et aid l2013bl : IMorandi fc Cuil 120141 ') : the gas den¬ 
sity slope reaches /3 ~ 0.9 at R200, in good agreement with 
the predictions of hydrodynamical simulations, but steeper 
than the values inferred from recent Suza ku observations 
l|Baut,z et al .11200?] : iKawaharada et aDl201Ct ). 


At present, Suzaku has mainly studied high mass clus¬ 
ters to the virial radius, primarily at low redshift, given 
its large PSF (~ 2'), which limits the size of the annuli 
used in the spatially resolved spectral analysis, reducing the 
angular resolution of the temperature and entropy profiles 
llBautz et alJ 120091 : lUrban et al.l l2013l l . It is clearly impor¬ 
tant to study clusters at higher redshift to investigate e.g. 
the evolution of the thermal properties of the ICM, and with 
independent X-ray observatories. 


In order to study the physical properties of the ICM in 
the virialization region we exploit the large archival dataset 
of Chandra clusters. We study the surface brightness pro¬ 
files of a sample of 320 clusters (0.056 < 2 < 1.24) with 


kT > 3 keV, observed with Chandra. We stackec0 the emis¬ 
sion measure profiles EM oc f n^dl of the cluster sam¬ 
ple to detect a signal out to and beyond ffeoo- We then 
measured the average emission measure, gas density and 
scatter in cluster outskirts. This pioneering work for the 
Chandra archive follows the original idea successfully ap 


Chandra archive toilows the original idea successfully ap¬ 
plied in ROSAT dEckert et al.112012 . see also|Dai_et_alj2003); 


and with the Sunyaev Zeldovich maps dPlagge et al 


20101 : 


IPlanck Collaboration et al .1 1201 .'l l, which allowed to detect 
the physical properties beyond 7?200- 

The paper is organized as follows. In [|2]we present the 
cluster sample, in m we describe the X-ray analysis and 
in m we discuss the systematics in the data analysis. The 
results on the gas density and gas mass fraction are pre¬ 
sented in m and m In Sj8] we summarize our conclusions. 
Throughout this work we assume the flat ACDM model, 
with matter density parameter = 0.3, cosmological con¬ 
stant density parameter £2 a = 0.7, and Hubble constant 
Hq = 100/ikm s _1 Mpc -1 where h = 0.7. Unless otherwise 
stated, we report the errors at the 68.3% confidence level. 


2 CLUSTER SAMPLE 


The sample we chose for this analysis is composed of 320 
clusters from the Chandra archive (Table [T]), with redshift 
range 2 = 0.056— 1.24 (the sample median redshift is ~ 0.4) 
and temperatures kT > 3 keV (the sample median temper¬ 
ature is ~ 7 keV). These sources have been observed with 
the ACIS-I imaging array. These clusters have been selected 
because the X-ray observations encompass R200 of each clus¬ 
ter, hence they are suitable for observations of cluster out¬ 
skirts. Moreover, source-free regions of the cluster observa¬ 
tion are always present, allowing measurements of the local 
background (Figure [U). The modeling of the background is 
indeed crucial for robust measurements of the physical pa¬ 
rameters in the outer regions. 

This sample encompasses most of the clusters 
from large sa mples studied i n pre vio us works, in¬ 
cludin g e.g. Vikhlinin et al] ( 20061'). Maughan et al.l 
J2008jh_ Cavagnolo et al.l ( 20091 4 . lEttori et al.l ( 2009a! ') and 
iMcDonald et al] (120131 ') . 

Although the sample was selected based on the quality 
of the existing observations and hence might be subject to 
selection effects, for the purpose of this work we did not 
require that the sample be complete, given the high level of 
self-similarity of the cluster outskirts. The X-ray properties 
of the galaxy clusters in our sample are presented in Table 
|T] while the number of net counts per bin for each cluster 
and the distribution of total net X-ray counts are shown in 
Figure IA1I 


1 We caution the reader that, from now on, with the word “stack¬ 
ing” we do not refer to an actual co-adding of the images, but 
rather to the median of the distribution of the radial emission 
measure profiles EM at radii R/R 200 , once re-scaled according 
the self-similar model (see [J3] for further details). 
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The galaxy cluster outskirts probed by Chandra 3 


Table 1. Table[l]is published in its entirety in the electronic edition of the Monthly Notices of the Royal Astronomical Society. A portion 
is shown here for guidance regarding its form and content. In the present table we report the X-ray properties of the galaxy clusters in 
our sample. For each object different columns report the cluster name, the redshift z, the exposure time t eX p, the identification number 
of the Chandra observations, the spectroscopic temperature T ew in the radial range 0.15 — 0.75 R 500 > the centroid shift (in units of R 500 ), 
and a flag for the presence or not of a cooling core (labeled CC and NCC, respectively). 


Cluster 

z 

^exp 

ks 

ID 

Tew 

(keV) 

w 

(xlO -2 ) 

CC/NCC 

A85 

0.056 

272.4 

904 4881 4882 4883 4884 4885 4886 4887 

4888 15173 16264 15174 16263 

6.45 ± 0.86 

0.49 ± 0.18 

CC 

A133 

0.057 

2393.0 

13442 13443 13444 13445 13446 13447 13448 13449 
13450 13451 13452 13453 13454 13455 13456 13457 
14333 14338 14343 14345 14346 14347 14354 13391 
13392 13518 3183 3710 12177 12178 12179 9897 

3.76 ±0.25 

0.26 ± 0.09 

CC 

A644 

0.070 

48.9 

2211 10420 10421 10422 10423 

8.12 ± 1.12 

1.56 ± 0.57 

CC 

A401 

0.075 

153.9 

10416 10417 10418 10419 14024 

5.44 ±0.86 

1.03 ± 0.38 

NCC 

A2029 

0.076 

29.2 

6101 10434 10435 10436 10437 

7.38 ±0.92 

0.21 ± 0.08 

CC 

A1650 

0.084 

218.4 

7691 5822 5823 6356 6357 6358 7242 10424 

10425 10426 10427 

5.89 ±0.86 

0.38 ±0.14 

CC 

A1068 

0.137 

19.5 

13595 13596 13597 13598 

6.38 ± 1.26 

0.29 ± 0.11 

CC 

A2276 

0.141 

39.5 

10411 

3.40 ±0.48 

0.49 ± 0.18 

NCC 

A1413 

0.143 

155.6 

5003 12194 12195 12196 13128 1661 5002 537 

7696 

7.78 ±0.22 

0.57 ±0.21 

CC 

A3402 

0.146 

19.8 

12267 

3.09 ±0.55 

1.39 ± 0.51 

NCC 


3 THE SELF-SIMILAR MODEL: MEASURING 

THE PHYSICAL PARAMETERS 

Stacking the emission measure of a cluster population is pos¬ 
sible thanks to the high level of self-similarity of the cluster 
outskirts. 

In this respect, the self-similar model (see, e.g.. iKaiser I 
Il98fil) gives a simple picture of the process of cluster for¬ 
mation in which the ICM physics is driven by the infall of 
cosmic baryons into the gravitational potential of the clus¬ 
ter DM halo. Power law scaling relations are expected under 
simplified models in which clusters are self-similar objects, 
having formed in single monolithic gravitational collapses 
and whose ICM is heated only by the shocks associated with 
the collapse. The underlying idea is that gravity is the only 
responsible for the observed values of the different physical 
properties of galaxy clusters; hence clusters are id entica l ob¬ 
jects when scaled by their mass or temperature (lLau et al.l 
120141) . 

Numerical simulations confirm that the DM compo¬ 
nent in clusters of galaxies, which represents the dominant 
fraction of the mass, has a remarkably self-similar behav¬ 
ior; however the baryonic component does not show the 
same level of self-similarity. For instance, deviation of the 
L — T relation in clusters with respect to the theoreti¬ 
cal value predicted by the previous scenario suggest that 
some energetic mechanis m, in addition to gravity, such as 
(pre)-heating and coolin g jBorg ani et al . 2005; B rvan II2000I : 
IMorandi fc Ettorl 120071 : Sunil 20 l3i ~ intervene to break the 
expected self-similarity of the ICM in the innermost re¬ 
gions. Consequently, the comparison of the self-similar scal¬ 
ing relations to observations allows us to evaluate the im¬ 
portance of the effects of non-gravitational processes on the 
ICM physics. While the central regions of clusters often ex¬ 


hibit complicated physical phenomena, such as AGN heating 
and cooling flows, the underlying physics of the outskirts is 
relatively simple, gravity-dominated and hence nearly self¬ 
similar. 

Assuming the spherical collapse model for the DM halo 
and the equation of hydrostatic equilibrium to describe the 
distribution of baryons into the DM potential well, in the 
self-similar model the cluster mass at an overdensity A (e.g. 
A = 500, 200) and temperature are related by: 

E z A 1 J 2 M tot oc T 3/2 ; (1) 

where the factor A z = Ax 

[1 + 82 (H 2 - 1) / (18tt 2 ) - 39 (n 2 - l) 2 / (18 tt 2 )] , with 
tl z = ± z) 3 /E 2 , acco unts for evol ution of clus¬ 

ters in an adiabatic scenario (IKaiser Hi). So we have 
Ra z oc (A//(p c , 2 A 2 )) 1/3 e x T^E-'A; 1 ' 2 . By setting 
f z = E Z (A Z /A) 1 ^ 2 (e.g. lEttori et al.l 120041 ). from the 
previous equations we can easily obtain the following 
relations: 

fzR.A z oc Tg/f 

/z-Mgas OC Tg/ S 2 (2) 

From the previous equation it follows that: 

EM oc J nidi oc M 2 a _ s /R 5 Az oc / 2 as f z T 3 ^ (3) 

assuming that all clusters have the same gas mass fraction 
/gas- Note that our sample is characterized by a relatively 
narrow temperature range (larger than 3 keV), such that 
the physical parameters of the clusters in our sample should 
show little dependence on gas fraction. The electron density 
n e (r) can be finally inferred by deprojecting the stacked X- 
ray emission measure profi le (Equat ion O via the spherical 
onion peeling method (IMorandi et al .112013al) . 
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Figure 1. Distribution of the clusters in our sample in a tempera¬ 
ture versus redshift diagram. The dashed line marks the region of 
the plot corresponding to 10 arcmin (the ACIS-I Chandra fov is 
16 x 16 arcmin). Note the some clusters, e.g. A133, A644, A2029, 
A1650 and A85, at low redshift 0.05), are characterized by val¬ 
ues of i ?200 larger than the fov of Chandra. These clusters have a 
large number of mosaic observations covering the outskirts. The 
blue and red dots correspond to clusters classified as cool core 
and non-cool core, respectively (see im . 


In our work we have considered all the physical quantity 
at fixed overdensity (A z = A), i.e. f z = E z in the above 
equations. 

In order to infer Raoo from the observables quantities, 
we rely on the existence of low-scatter scaling relations be¬ 
tween the total thermal energy Yx = M gaa Tx and the to¬ 
tal mass, as predicted by self-similar theory and confirmed 
by high-resolution cosmological simulations (iKravtsov et al.l 
2006). In this respect we used the Yx —M relation measured 
for the lSun et al.1 (120091) sample of clusters: 


E 2/5 M 5 oo 



5 B i 
2 ± 

1.63 x 10 14 


Y x 

3.60 x !O 13 M 0 keV 


0.571 

(4) 


with B = 0.571. Hence we recovered Rsoo- 

We point out that we determined Tx in the radial range 
0.15 — 0.75 J? 5 oo 1 H4.2I) . while the previous relation has been 
calibrated in the range 0.15 — 1 Rsoo- Our choice is moti¬ 
vated by the fact that this is a straightforward measurement 
for our intermediate and high-z clusters because exposures 
were designed to provide a sufficient statistical accuracy to 
infer the spectral temperature in the aforementioned radial 
range. We correct our Tx measurements in order to match 
the definition employed in Equation [4] via mock simulations 
by assuming our average gas density profile (Sj 6 j and the 
averag e temperature profile as measured bv lVikhlinin et ahl 
|200g). This correction is 0.96 ± 0.02. 

i ?500 and Yx were computed iteratively, measuring the 
temperature in the aperture 0.15 — 0.75 R 500 and the gas 
mass within R 500 , computing a new Yx, and hence estimat¬ 
ing a new value of R 500 . 


In order to infer Ra (A = 500, 200,100) from the known 
value of -R 500 , we apply the following statistical method. 
First, we assumed a NFW distribution with parameters 
(c, i? 20 o); second, we inferred a value^of the concentrat ion pa¬ 
rameter via the c — M v i r relation of lDuffv et al.1 (l2008li based 
on the results from numerical simulations. This provides an 
estimate of Ra for given value of -R 500 . For our sample we 
have R 500 = (0.65±0.01)J?2oo and R 100 = (1.36±0.01)7?2oo 
The gas fraction within an overdensity A can be com¬ 
puted directly from the gas density profiles: 

/gas,A = Mg J’ A = ~T~ [ Pgas(x)x 2 dx. (5) 

M a Ap c , z J 0 

with x = R/Ra- 


4 X-RAY ANALYSIS 

D escription of the X-ray a nalysis methodology can be found 
in IMorandi et ali (l2013bll . Here we briefly summarize the 
most relevant aspects of our data reduction and analysis. 


4.1 X-ray data reduction 

All data were reprocessed from the level 1 event Hies us¬ 
ing the Cl AO data analysis package - version 4.6.1 - and 
the calibration database (CALDB 4.6.3) distributed by the 
Chandra X-ray Observatory Center. Note that the ACIS QE 
contamination model, which has been released in CALDB 
and is associated with the deposition of one or more mate¬ 
rials on the ACIS detectors or optical blocking filters, is the 
latest available (version N0009). 

We analyzed 613 datasets, corresponding to 320 galaxy 
clusters, retrieved from the NASA HEASARC archive with 
a total exposure time ~ 20 Ms for the whole sample. All the 
observations are carried out by using the ACIS-I CCD. 

We reprocessed the level-1 event files to include the ap¬ 
propriate gain maps and calibration products. As part of the 
data reduction, corrections were made for afterglows, charge 
transfer inefficiency (CTI), bad pixels and solar flares. We 
used the acis_process_events tool to check for the pres¬ 
ence of cosmic-ray background events, correct for spatial 
gain variations due to charge transfer inefficiency and re¬ 
compute the event grades. Then we filtered the data to in¬ 
clude the standard events grades 0, 2, 3, 4 and 6 only, and 
therefore we filtered for the Good Time Intervals (GTIs) 
supplied, which are contained in the fIt 1. fits file. 

A careful screening of the background light curve is nec¬ 
essary to discard contaminating flare events. A common way 
of removing these periods of particle background flares is 
to create a lightcurve of the local background, once point 
sources of high and variable emission are excluded. The 
lightcurve is then created by using the tool dmextract. In 
order to clean the datasets of periods of anomalous back¬ 
ground rates, we used the deflare script, so as to filter out 
the times where the background count rate exceed ±3er of 
the mean value. Most observations were taken in VFAINT 
mode, and in this case we applied VFAINT cleaning to both 
the cluster and blank-sky observations. Finally, we filtered 
the ACIS event files on energy selecting the range 0.3-12 
keV, so as to obtain a level-2 event file. 

Point sources and extended substructures were detected 
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and removed using the script wavedetect, which provides 
candidate point sources, and the result was then checked 
through visual inspection. 

4.2 X-ray surface brightness analysis 

We produced X-ray images from the level-2 event file. Our 
goal is indeed to measure the emission measure and the gas 
density profile in a non-parametric way from the surface 
brightness. Multiple observations were reduced individually 
to apply the correct calibration to each dataset. The cluster 
surface brightness profile is obtained from merged images. 

The X-ray images were extracted from the level-2 event 
files in the energy range 0.7 — 2.0 keV. This energy band was 
chosen in order to: i) maximize the S/N in the outskirts of 
the clusters, and ii) minimize the dependency of the cool¬ 
ing function on the temperature and metallicity. Concern¬ 
ing the latter point, the integrated cooling function in the 
(0.7 — 2) keV band is approximately given by A (T) oc T~ a , 
with —0.02 < a < 0.2 for T 2 keV (in all the clusters the 
temperature is not expected to fall below ~2 keV), such that 
systematic uncertainties in the estimated projected temper¬ 
ature have a negligible impact on the gas density. In partic¬ 
ular, we point out that we inferred the density profiles out 
to -Rioo from the stacked emission measure profiles. 

We then corrected the images by the exposure maps 
to remove the vignetting effects. We created an exposure- 
corrected image from a set of observations using the 
merge_obs tool to combine the ACIS-I observations. All 
maps were checked by visual inspection at each stage of the 
process. 

Next, we determined the centroid ( x c , y c ) of the surface 
brightness image by locating the position where the deriva¬ 
tives of the surface brightness variation along two orthogonal 
(e.g., X and Y) directions become zero. 

Since for all the clusters the field-of-view encompasses 
the X-ray boundary (i?ioo), the emission does not extend 
across the entire detector, and we have a sufficiently large 
region to obtain a local background. In all the clusters the lo¬ 
cal background was measured from a region beyond 1.5 R 200 , 
where the surface brightness was approximately constant, 
i.e. it has reached the background level. 

The local background is a combination of a particle 
component that is not vignetted, and a sky component that 
is vignetted. To determine the surface brightness of the clus¬ 
ter and of the local soft X-ray background, an accurate 
procedure consists of subtracting the non-vignetted parti¬ 
cle component as measured from Chandra observations in 
which the ACIS detector was stowed, after rescaling the 
stowed background to match the 9.5 — 12 keV cluster count 
rate of the observations (where the Chandra effective area is 
negligible and thus very little cluster emission is expected). 
Indeed, although the partic le back ground flux may vary with 
time, iHickox fc Markevitchl (l2006li has shown that the spec¬ 
tral distribution of the particle background is remarkably 
stable, even in the presence of changes in the overall flux, 
and that the ratio of soft-to-hard (2-7 keV to 9.5-12 keV) 
count rates remains constant to within < 2 percent in time. 
We verified that all our sources are characterized by remark¬ 
ably stable background by comparing the spectra of the local 
background and of the blank sky fields. 

In order to correct for vignetting, a common approach 


is to divide the counts image by an exposure map, weighted 
according to a specific model for the incident spectrum, to 
rescale all parts of the image to the same relative exposure. 
Nevertheless, since the exposure map is both energy and 
position dependent, we cannot assume the same exposure 
map for source and background emission. Our model of the 
desired source surface brightness Sx,s breaks down the ob¬ 
served emission Sx,obs into four components: 

5x, 0 i, s (x) oc Sx,s(x.) ■ Cs(x, E) + Sx,bkg(x) ■ Cbkg{x,E) 

+Sx, P b{x, E(t)) + SxA X, E) (6) 

where Sx,s{x) is the source surface brightness in the pixel 
x = ( x,y ), Sx,bkg is the cosmic X-ray background (CXB) 
background and Sx, P b is the (time dependent) particle back¬ 
ground. C„ (x, E) and Cbk g (x, E) refer to the exposure maps 
for source and background, respectively. 

Sx,r(x, E) refers to the readout artifact background. 
We simulated, re-normalized and accounted for this read- 
out backg round usi ng the make_readout_bg routine (see 
IHickox fc Markevitchl ( 20061 ) for further details). 

The radial surface brightness profile is thus derived with 
the exposure correction and particle background subtraction 
using the scaled stowed background. The region where the 
CXB is more dominant than the cluster emission can be de¬ 
termined from the flattened portion at the outer region of 
the surface brightness profile. We applied a direct subtrac¬ 
tion of the CXB+particle background by means of Equation 
[6] In order to calculate Cbk g {x, E), we modeled the soft CXB 
component by an absorbed power law with index 1.4 and 
two thermal components at zero redshift, one unabsorbed 
component with a fixed temperature of 0.1 keV and another 
absorbed component with a temperat ure derived from spec¬ 
tral fits (~ 0.25 keV. ISun et al.ll2009h . 

The boundary radius of the X-ray and surface bright¬ 
ness analysis corresponds to a ratio of source to background 
flux ~ 15-40 percent at i? 200 - Thus, a careful analysis of the 
systematics and how these uncertainties propagate into the 
determination of the physical parameters is required ( 1 )5.21) . 

With robust background subtraction and modeling, 
from the X-ray images we measure the emission measure 
profile EM oc f nidi and calculated the median distribution 
to compute stacked profiles. We then scaled the emission 
measure profiles and the deprojected density profiles for the 
clusters in our sample. A self-similar scaling was applied to 
the emission-measure profiles, i.e. each profile was rescaled 
by the quantity El(kT/ 10 keV) 1 / 2 (§3]), following their evo¬ 
lution with redshift and dependency on the cluster mass, 
since the profiles is expected to show a remarkable level of 
self-similarity outside of the core (R > 0.2J?20o)- From the 
deprojection of the median emission measure profile we then 
recover the expected density profile. 


4.3 X-ray spectral analysis 

The global temperatures have been recovered by fitting spec¬ 
tra from different observations simultaneously. The spec¬ 
tral analysis was performed by extracting the source spec¬ 
tra from circular annuli around the centroid of the surface 
brightness and by using the Cl AO specextract tool from 
each observation. The spectral fit was performed by simul¬ 
taneously fitting an absorbed thermal emission model in the 
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energy range 0.7-7 keV and in th e reg ion 0.15 — 0.75 -R500. 
We used the XSPEC package (lArnau Jl996l . version 12.8.2) 
to perform th e spectral fi t. We adopted the APEC emis- 
sivity model llFoster et al]|20l5 ) and the AtomDB (version 
2 . 0 . 2 ) database of atomic data, and we employ ed the so¬ 
lar abundance ratios from lAsplund et al.l (200a) • We also 
used the Tuebingen-Boulder absorption model (tbabs) for 
X-ray absorption by the ISM. Free parameters in the APEC 
model are temperature, metal abundance, normalization. 
We fixed the hydrogen column density Nh to the Galactic 
value by using the Lei den/Argentine/Bonn (LAB) Hl-survey 
dKalberla et al.l 120051') . The redshift to the value obtained 
from optical spectroscopy. We also group photons into bins 
of at least 20 counts per energy channel and applying the 
X 2 -statistics. We analyzed the observations individually to 
check for consistency before jointly fitting multiple datasets 
referring to a single cluster. The background spectra have 
been extracted from regions of the same exposure for the 
ACIS-I observations, for which we always have some areas 
(> i?ioo) free from source emission. 

We also checked for systematic errors due to possible 
source contamination of the background regions, which has 
been estimated by considering the predictions of the sur¬ 
face brightness from hydrodynamical numerical simulations 
including cooling^ star formation and supernovae feedback 
(lRoncarelli et al]i200a j. The contamination is always of the 
order a few percent, and this translates into a negligible bias 
of the determination of the spectral global temperature. 

We finally validated the method to recover the global 
temperature for each cluster by creating a mock spectrum 
via a single-temperature absorbed thermal model with pa¬ 
rameters fixed to those from the aforementioned spectral 
analysis, and then adding a background spectrum from a 
region free of the source emission. We thus repeated the 
spectral analysis on the mock datasets and we found a very 
good agreement between input and recovered spectral pa¬ 
rameters. 

I 11 the appendix we present a comparison between 
spectral temperatures and metallicities recovered via the 
AtomDB databases 2.0.2 and 1.3.1, via different models of 
X-ray absorption by the ISM, and by accounting also for 
the molecular and ionized hydrogen in the hydrogen column 
densities. 

4.4 Morphological type classification 

In this section we describe how the clusters were classified 
based on their morphological type, i.e. based on the presence 
of a cool core and/or dynamical state of clusters. Dividing 
the data into these subsamples allows us to investigate the 
effect of cool cores and morphological disturbance on a clus¬ 
ters position with respect to the mean relation. 

We adopted the centroid shift method as our method 
to in fer to the dynamical state of clusters dMaughan et al.l 
I 2008 T ). Centroid shifts have been measured from the vari¬ 
ation of the centroid of cluster emission within apertures 
of increasing radii Ri, which have similarly been previously 
used to infer the emission measure. The centroid shift, w, 
was defined as the standard deviation of these centroids in 
units of -R500. 

Centroid shifts were measured from exposure-corrected 
images in order to eliminate the effects of vignetting. We 
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Figure 2. Plot of centroid shifts versus electron density at 
0.03 -R 500 for the whole sample. 


masked out point sources, chip gaps and all the substruc¬ 
tures which have been similarly excluded in recovering the 
emission measure profile. This generates a 2D mask I(x — 
x c ,y — yc), with x c ,yc being the centroid as determined in 
1 14.21 We then symmetrized this mask, i.e. I(x — x c , y — yc) = 
I(—x + x c, — V + Vc), in order to remove any spurious depen¬ 
dency of the centroid shifts on the adopted masking, i.e. a 
perfect azimuthally symmetric source must generate a zero 
centroid shift regardless of the adopted mask. This definition 
adds the benefits to robustly relate w to the global dynam¬ 
ical state of clusters, while eliminating the impact of point 
sources, chip gaps and substructures similarity masked in 
the emission measure profiles. 

The measured centroid shifts are summarized in Table 
□ The errors quoted for the w are the statistical uncer¬ 
tainties on a standard deviation calculated from n measure¬ 
ments (w^/2/(n — 1)). Objects with center shift parameter 
w > O.OIR 500 are classified as morphologically disturbed. 

Next, systems are classified as cool cores (CC) or 
non-cool cores (NCC) on the basis of density E~ 2 n e ,0 at 
0.03 7?5 oo- The threshold we use to define a cool core systems 
is E~ 2 n e fi > 1-5 x 10~ 2 cm~ 3 . We compa red our definition of 
CC systems with lCavagnolo et al.l (l2009l l in m whose sam¬ 
ple shares 93 objects in common with ours. They use the 
entropy threshold Kq = 30 — 50 keV cm 2 to approximately 
demarcate the division between CC and NCC. 

Note that centroid shifts and central density are broadly 
anti-correlated (Figure[2]), being both a proxy of the dynami¬ 
cal state of a cluster, i.e. cool-core system have larger central 
density and overall more relaxed dynamical state. 

We finally measured the flattening and orientation of 
the X-ray surface brightness. We computed the moments 
of the surface brightness within a circular region of ra- 
dius -R 500 cent e red on the centroid of the X-ray image (see 
iMorandi et al.l (I2010T ) for further details on this method). 
This allows us to estimate the ICM eccentricity on the plane 
of the sky and the orientation (position angle) of an elliptical 
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X-ray surface brightness distribution. We use this position 
angle as a proxy of the direction of the large-scale filament 
(see discussion in 


5 ANALYSIS OF THE SYSTEMATICS IN THE 
DATA ANALYSIS 

5.1 Validation of the stacking method 

We validated our stacking method via simulations of mock 
datasets including systematics. From the rescaled emission 
measure, averaged on the whole sample, we created mock 
azimuthally-averaged source brightness images Sx,s(r) for 
each observation, according to the observation aspect in¬ 
formation. We then used blank-sky fields to produce realis¬ 
tic backgrounds S^ 9 (x) for each observation. We extracted 
background images from the blank-field background data 
provided by the ACIS calibration team in the same chip 
regions as in the observed cluster. The blank-sky observa¬ 
tions event files underwent a reduction procedure consistent 
with that applied to the cluster data, after being reprojected 
onto the sky according to the observation aspect informa¬ 
tion by using the reproject_events tool. We then added 
these background images to the source brightness. In or¬ 
der to estimate the contribution of the particle background, 
we used stowed background observations, renormalized to 
match the blank-sky count rate at high energy (9.5-12 keV). 
We point out that in this way we capture biases due both 
possible spatial variation of the CXB and temporal varia¬ 
tions of the particle background spectrum. We finally ap¬ 
plied the whole procedure described in H4.ll fEauation[51) on 
the mock observations. A comparison between the true and 
measured emission measure profile reveals that the distri¬ 
bution of the residual (normalized by the measurement un¬ 
certainties) follows the expected Gaussian distribution with 
zero mean, with reduced \ 2 of the order of the unity. This 
indicates that the measurement errors fairly represent the 
error budget and our recovered physical parameters are not 
significantly biased. 


5.2 Systematics in background modelling 

Accurate measurements of the ICM in the cluster outskirts 
require an accurate knowledge of the components of the X- 
ray background, and an accurate understanding of how the 
uncertainties of the background modeling propagate on the 
desired physical parameters. 

For more details and a thorough validation of the 
method to g auge statistical and sy stematic errors, we refer 
the reader to lMorandi fe Cuil (l2014h . Here we briefly outline 
the key aspects. 

The cosmic X-ray background (CXB) consists of the lo¬ 
cal hot bubble emission, the thermal emission of the galac¬ 
tic halo and the contribution of unresolved point sources 
(mostly AGNs). The latter component is modeled as a power 
law with index 1.4. When analyzing regions of finite size 
we must account for Poisson fluctuations of these unre¬ 
solved point sources. The expected deviation from the av¬ 
erage value for a given observed solid angle S2 resolved to a 


threshold flux SWes is: 

fS thres / J\J \ 

<XCXB = (I/O) J (^) x S 2 dS (7) 

where dN/dS is the c umulative flux distribution of point 
sources as employed bv lMoretti et al.1 (|2003h . The threshold 
flux Sthres has been calculated by determination of the local 
flux limit to which we can robustly identify a point source, 
com monly known as sensitivity map (see, e.g.. lEhlert et all 
120131 . for further details on the method). Our joint Chandra 
observations of the outskirts allow the CXB to be resolved 


soft band. Solving Equation [7] this translates into an error 
in the range <)3% on the surface brightness in the outer 
volumes for the average emission measure, while it becomes 
negligible in the inner volumes. If we consider the subsample 
of clusters at higher redshift (z > 0.3), this error becomes 
appreciable (< 10 percent in the outer volumes). 

We point out that the previous analysis does not take 
into account the possible systematics due to variation of the 
PSF across the ACIS-I detector, the PSF being about one 
order to magnitude larger (~ 10 — 20 arcsec) at the edges 
of the CCD with respect to the aimpoints (~ 1 arcsec); 
neither it accounts for the different contamination of the 
galaxy cluster signal in detecting point sources across the 
CCD. Indeed, the signal of point sources near the edge of 
the field would be diluted on larger angular scales due to 
the larger Chandra PSF, but at the same time would be less 
contaminated by the cluster signal, the cluster center be¬ 
ing observed customarily near the center of the CCD. The 
first (second) effect would make more difficult (easier) to de¬ 
tect point sources with respect to the CCD center via the 
wavelet decomposition algorithm described in H4.ll In order 
to gauge the impact of these systematics we perform Monte 
Ca rlo (MC) sim u lation s of the point sources as employed 
by iMoretti et al.l (l2003l l , which should account for most of 
the CXB background. We then added particle and galactic 
background, which has been convolved with the Chandra 
PSF at the position (x, y) of the CCD. We finally added 
mock cluster signal (see SJ3] for details on the mocks). For 
each MC simulation we extracted point source by mimick¬ 
ing the procedure of wavelet detection algorithm. This allow 
us to estimate the bias (downwards) due to the aforemen¬ 
tioned systematics, which is < 3% (4%) on the gas density 
at R, 20 o{Rioo)- 

We then estimated further systematics in the data anal¬ 
ysis due to background subtraction. For this purpose, we use 
stowed background and local background taken from source- 
free regions of the cluster observations in order to model 
the contribution of the CXB and particle background, re¬ 
spectively. The two major sources of systematic uncertain¬ 
ties are: i) the renormalizion of the stowed background to 
match the cluster count rate at high energy (9.5-12 keV) 
due to Poisson errors; and ii) uncertainties in the parti¬ 
cle background across its spectrum. Note that the prop¬ 
erties of the particle background are well known via the 
stowed background observations, with uncertainties in the 
back ground across its spect rum which are within £ 2 per¬ 
cent (IHickox fc Markevitcbl (l2006l 'l and see also discussion in 
ISun et al. ( 2009T) 1. 


Finally, we estimated the uncertainty in our measure¬ 
ments due to soft X-ray background variations across the 
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field of view. We divided the regions used to obtain the 
background model into independent subregions with differ¬ 
ent azimuthal direction with respect to the cluster center. 
We implemented a bootstrap approach, where we randomly 
picked these background subregions, allowing repetitions, 
and re-determined the background level via Equation [6] We 
repeated this procedure 10 s times, and determined the dis¬ 
tribution of the average background. The uncertainties due 
to CXB fluctuations are always smaller than the statisti¬ 
cal error bars determined by assuming that the soft X-ray 
background is spatially constant. 

All the aforementioned uncertainties on the background 
modelling (cosmic X-ray background, particle background, 
soft X-ray background variations and unresolved point 
sources) have been propagated in recovering the desired 
physical parameters, by means of MC randomizations of the 
background including both statistical and systematic uncer¬ 
tainties. 

5.3 The impact of undetected clumps, gas 
inhomogeneities and multi-temperature 
distribution 

The CDM scenario predicts a picture where clumps at sub- 
virial temperature are infalling along filaments and accret¬ 
ing onto the cluster outskirts. Given the low signal-to-noise, 
low sensitivity/integration time and/or poor spatial resolu¬ 
tion of the X-ray telescopes, it might be impossible to mask 
out cold substructures e.g. in the X-ray analysis, and they 
might remain undetected. A cold phase (T < 1 keV) of the 
ICM would also coexist in the outskirts with the hot ICM 
at virial temperatures (> 2 — 4 keV), leading to a multi¬ 
phase structure of the ICM. The presence of a multiphase 
structure of the ICM will bias the spectroscopic tempera¬ 
tures and gas density. Spectral temperature measurements 
indeed customarily hinge on a single-temperature absorbed 
thermal model, since the data quality is in general not suf¬ 
ficient to constrain a multitemperature model, especially in 
the outskirts where the signal is overwhelmed by the noise. 

The presence of undetected dense and cold clumps at 
temperatures 0.5 <J kT <) 2 keV in the ICM will typi¬ 
cally enhance the X-ray emission, since the instrument is 
very sensitive to the low energy features of these clumps 
whose exponential bremsstrahlung cutoff will be at energies 
~ kT/{ 1 + z). Therefore, besides upon the dynamical state 
of the cluster, the impact of dense clumps will also depend 
on the redshift of the observation, the Chandra effective area 
being small for photon energies smaller than 0.7 keV. Thus, 
if this dense clumps at subvirial temperatures are very cold 
(< 1 keV), high-z clusters might be less biased by the pres¬ 
ence of a gas at low temperatures, but at the same time it 
will be more difficult to detect these substructures given the 
cosmological dimming. Shallow observations are also more 
prone to be biased by undetected clumps. In Figure [3] we 
present the bias on the gas density due to undetected sub¬ 
structures at subvirial temperatures. 

The qualitative picture emerging by this toy model is 
that the presence of undetected substructures might bias 
upwards the measured gas density, with a magnitude larger 
for shallow observations and low redshifts. Nevertheless, by 
dividing our sample into two subsamples, clusters with to¬ 
tal net counts smaller and larger than the median value of 



Figure 3. Bias (upwards) on the gas density at -R 500 due to un¬ 
detected substructures at subvirial temperatures T~ 0.4 keV and 
0.6 keV. The plotted redshift corresponds to the range where the 
subclumps fall below the detection threshold of the observation 
(JOJ, i.e. they are undetected. We assumed a background level 
of 2 X 10 —11 erg cm -2 s~ 1 deg’ 2 in the soft band, a filling factor 
of 1% for these subclumps, a cluster temperature of 4 keV and 
self-similar gas density profile for both cluster and subclumps. 
The solid (dashed) line refers to an observation time of 100 (10) 
ks. 


the whole sample, respectively, and repeating our stacking 
procedure, we do not find any significant difference in the 
recovered physical parameters, which might arise given the 
different impact of substructures upon the observation time. 
While the superb angular resolution of Chandra is essential 
for detecting small-scale clumps and distinguishing the X- 
ray emissions arising from clumps and diffuse components, 
we caution the reader that it possible we might actually 
missing some substructures even in deep observations, given 
their low signal-to-noise in the outskirts, which could en¬ 
hance the X-ray emission and hence the gas density. 

Moreover, inhomogeneities in the gas distribution can 
also lead to the overestimate of the observed gas density 

llNagai fe LaullioiH: IZhuravleva et al.|l2013l: iRoncarelli et al] 

120131) , which in turn introduces biases in global cluster prop¬ 
erties, such as the gas mass fraction (IBattaglia et al.ll201fil ) 
and the low-scatter X-r ay mass proxy, such as Yx = M gas Tx 
llKhedekar et al] 120131 ). At R 200 hydrodynamical simula- 
tions predict a gas clumping factor in the range ~ 1.3 — 2 
llNagai fe haul lioTil : IZhuravleva et al.|l2013l : IBattaglia et alJ 
120131) . Gas clumping factor inferred from th ese sim ulati ons 
is in reasonable agreement wi th observations llMorandi et al.1 
l2013bl : iMorandi fe Cuill2014l ). If they are not properly un¬ 
derstood and modeled, these non-equilibrium processes can 
limit the use of galaxy clusters as cosmological probes. 

Finally, in order to understand how a multiphase struc¬ 
ture of the ICM can bias X-ray observations, we present 
a simple toy problem and we analyze the effect of multi¬ 
temperature distribution on the physical observables (i.e. 
X-ray spectroscopic temperature and gas density). A cold 
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Figure 4. Bias on the gas density due to a multi-temperature gas distribution of the ICM. We fitted an absorbed thermal emission 
model at a single emission temperature T spe c on mock spectra of a two-phase plasma at temperatures 7) Kli and T co i<j and in pressure 
equilibrium, with gas mass fraction of the cold gas component //. We assume that the gas mass-weighted temperature of this two-phase 
plasma (as measured e.g. by SZ) is constant, i.e. T mw = 4 keV. The left (right) panel refers to an observation redshift of 0.05 (1). 


phase of the ICM would systematically bias down X-ray 
spectroscopic temperatures with respect to the average (gas 
mass-weighted) temperature e.g. probed through SZ. We 
consider a two-phase gas at temperatures T co id and Thot, 
with gas masses m co id and mhot, respectively, and in pres¬ 
sure equilibrium. Hence the gas mass-weighted temperature 
T mw reads: 

T tnw — (1 /t ) / Ik,I A//'/cold , jl = Ttl cold/(tUcold“fUthot) (8) 

with T mw = T mw (Thot, Tcoid, m) i T co id < T mw < Thot- 

We generated mock spectra of the previous plasma by 
means of absorbed APEC emission models at temperatures 
Thot and T co id, and under the boundary condition that gas 
mass-weighted temperature (Equation [8| is constant, i.e. 
T mw = 4 keV. The hydrogen column density Nh is fixed 
to 10 20 cm -2 . The current emission model is then folded 
through response curves (ARF and RMF) of Chandra. 

Since T mw is constant, the problem to recover the spec¬ 
tral temperature T spec of the two-phase plasma is fully de¬ 
scribed via two independent variables, e.g. T co id and /r, while 
Thot is a dependent variable which can be inferred via Equa¬ 
tion [8j Thus, T sp ec = Ts P ec(Tcoid, ft), where we measured the 
spectral temperature by fitting an absorbed APEC emis¬ 
sion models at a single-temperature T ape c on the mock spec¬ 
tra of the two-phase plasma. The final goal is to infer the 
bias on the gas density, which arises since the spectroscopic 
temperature T spe c is then used to convert observed X-ray 
counts into an emission measure (and thus gas density) via 
a single-temperature absorbed thermal model. In Figure 2] 
we present the results of the bias on the gas density for a 
grid of values of T co id and fr. 


We can observe that the gas density on the two-phase 
plasma is biased with respect to the true density. Since the 
primary emission mechanism in X-ray clusters is thermal 
bremsstrahlung and emission lines, the spectral fit is indeed 
guided, primarily, by the exponential bremsstrahlung cut¬ 
off at high energies, the iron line at energies ~ 7/(1 + z) 
keV, and the lower energy (<) 2 keV) emission lines. Com¬ 
plications arise since the true spectrum of the astrophysical 
source must be convolved with the instrumental response of 
the X-ray telescope. Chandra effective area, for example, is 
small for photon energies smaller than 0.6 keV or greater 
than 7 keV. Hence, the presence of a cool phase at temper¬ 
ature 0.5 <) T <) 2 keV in the ICM will typically bias down 
the spectroscopic temperature (single-temperature absorbed 
thermal model), since the instrument is very sensitive to the 
low energy features, such as low energy emission lines. In 
Figure 2] we can see that the gas density is biased upwards 
as long as ~ T co id/(1 + z) )> 0.25 keV, a with non-negligible 
dependency on the observation redshift, fr and T co id- Coun¬ 
terintuitively, for low temperatures ~ T co id/(l + z) < 0.25 
keV and larger values of the cold phase fraction fi )> 0.2, the 
measured gas density is biased downwards of ~ 0—20%. This 
is due to the fact that most of the emission captured by the 
X-ray telescope arises from the hot phase of the ICM, with 
temperature Thot )> T mw for T co id/(l + z) <, 1 keV. On the 
contrary, most the emission of the cold phase would be at 
very low energies, where the Chandra effective area is small, 
and related to emission lines rather than bremsstrahlung. 
This bias is increasing for higher observation redshifts. 

A multiphase structure of the ICM would retain signa¬ 
tures of the ’melting pot’ in the virialization region, where 
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infalling clumps of matter and gas are predicted to have, 
typically, higher density and cooler temperature than their 
surroundings. Current measurements via X-ray telescopes, 
which customarily hinge on a single-temperature absorbed 
thermal model, can be thus biased. Future instruments (e.g., 
the upcoming Astro-H and an envisaged Athena-like mis¬ 
sion) will have enough collecting area and spectral resolu¬ 
tion to better distinguish multi-temperature structures spec¬ 
trally. Moreover, more detailed simulations (e.g. hydrody- 
namical simulations) including these effects might be needed 
to address the impact of these biases. 

5.4 Other sources of systematics 

We checked the effect of the presence of a temperature pro¬ 
file on the measured emission measures. The conversion from 
ACIS-I count rate to emission measure has been determined 
via an absorbed thermal emission model folded with the 
AICS-I response. This conversion factor was inferred for each 
cluster by using the global spectroscopic temperature in the 
range 0.15 — 0.75 -R 200 rather than the temperature profile. 
As we previously 1 >14.211 pointed out, the conversion from 
ACIS-I count rate to emission measure is highly insensi¬ 
tive to the temperature (for T > 2 keV). The conversion 
factor changes at most by 8% for low-temperature (T ~ 3 
keV) clusters, which translates into a bias 4% on the gas 
density, while it is negligible (<j 1%) for hot (T )> 6 keV) 
systems. To further check the effect of the presence of a 
temperature profile on the measured stacked emissivity and 
hence gas density, we generated mock event files for each 
cluster under these assumptions: i) the gas density distribu¬ 
tion for each cluster is self-similar and fixed to that mea¬ 
sured via our stacking method as described in > 14. 21 ii) we 
assume the expression for the self-similar temperature pro¬ 
file that lVikhlinin et al.l (120061 1 showed to reproduce well the 
temperature gradients in nearby relaxed systems, with nor¬ 
malization fixed to the spectral temperature of each source. 
We repeated the stacking procedure described in > 14.2l bv as¬ 
suming first isothermality and then by accounting for the 
temperature profile. We verified that the recovered physical 
parameters are biased downwards in a negligibly way (<J 2% 
for the emission measure) by our method. 

As a further test, we compared our gas mass at R 500 
(recovered via deprojection of individual em i ssion measure 
profiles) with the results of lMaughan et ahl (I2006T ). whose 
sample sha res 9 5 o bjects in common with ours. We also com¬ 
pared with lMartino et al.l d2014l l. whose sample contains 32 
objects in common with us. We found no significant differ¬ 
ences between the recovered gas masses (Appendix [D]). We 
point out that these individual gas masses profiles have been 
recovered only for the purpose of this comparison, since we 
actually use stacked profiles across the paper. 

Next, we compared the stacked emission measure pro¬ 
file by using median and weighted average of the emission 
measure profiles, and we did not find any statistically signif¬ 
icant difference for R :> O.Q 0 R 200 . Similar conclusions hold 
for the gas density (Figure 0- This is a consequence of the 
fact that the errors (measurement and systematic) are sym¬ 
metric. However, a weighted-average would significantly un¬ 
derestimate the final uncertainties on the stacked emission 
measure, since it accounts only for the measurement errors 
but it neglects the intrinsic scatter of the emission measure 


profiles. Hereafter, we will refer to the results recovered via 
the median of the emission measure profiles. 

Finally, we estimated the b ias due to uncertainty in 
the M 500 — Fx.soo relation t|Sun et alil2009) i on the desired 
physical parameters due to the presence of non-thermal 
pressure. The latter can indeed bias downwards the hy¬ 
drostatic masses through which the M500 — Yx, 500 relation 
has been calibrated. We define the mass bias b between 
the true and hydrostatic mass Mhydro , with both masses 
defined at a fixed density contrast of 500. Hence we have 
that Mtrue = Mhyd.ro /(1 — 6). We further assume an univer¬ 
sal gas density profile as determined in our work (>0, and 
the expression f or the self-similar temperature profile from 
IVikhlinin et al.l d2006l !. Thus, by changing the bias param¬ 
eter b and correcting accordingly the normalization of the 
A /500 — Yx ,500 relation, we can infer the R 500 = R 50 o(T,b) 
relation from Equation [4] which can adequately described 
by the following expression: 

eSv) M “" k - < 9 > 

The bias on the true gas density becomes n e ,true/n e ,mens — 
(J?50o(6)/7?500,meas) _ “. Note that the bias parameter b 
refers to the non-thermal pressure support at R500, since 
the M500 — Va',5oo relation has been calibrated at this over¬ 
densities. For larger radii we do not expect a larger impact 
of the non-thermal pressure, since we used the Yy.soo es¬ 
timator in order to infer R500; hence R200 and R100 have 
been statistically determined by assuming an NFW distri¬ 
bution for the total matter (>0. Fro m Equation 0 and for 
b ~ Q l (see, e.g. i Mahdav i et ahlfioOSl : iMorandi fc Limousinl 
l2012l:lMorandi et al.ll201^ the presence of non-thermal pres¬ 
sure can bias upwards the gas density of ~ 15(25)% at 
7?50 o(7?2oo). 


6 THE OUTSKIRTS OF CLUSTER: GAS 

DENSITY 

Results for the gas density, gas density slope /3 = 

— l/3dlog(n e )/dlog(r) and gas fraction / gas (< R) = 
A/gas(< R)/Mtot{< R) for our sample are presented in Ta¬ 
ble [2] We divided the data into subsamples to investigate 
the dependency on the redshift z, cluster temperature T, 
cluster morphological type (relaxed/unrelaxed). We conven¬ 
tionally defined relaxed (unrelaxed) system with w < 10~ 2 
(w > 10 -2 ). We remember that the centroid shifts and 
central density are anti-correlated (Figure 0, being both a 
proxy of the presence of a cool-core: thus the chosen thresh¬ 
old of w roughly demarcates the division between CC/NCC. 

The quoted uncertainties reflect the total error budget 
(measurement error + intrinsic scatter) by means of MC 
randomizations. 

In order to gauge the bias implicit in our spherical mod¬ 
eling, we further repeated the X-ray analysis by restricting 
the stacking procedure ( > 14.211 in two sectors (/1 and f 2 in 
Table 0. The sector /2 has position angle fixed to the ma¬ 
jor axis of the X-ray brightness distribution and with an 
aperture of 60 degrees, i.e. encompassing the cosmic fila¬ 
ment (see > 14.41 for further discussion), while the other (/ 1 ) 
characterizes the remaining volumes of the cluster. 


R 500 = E - 1 1095 
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Figure 5. Left: The Chandra emission measure profiles of all clusters in the sample, re-scaled according the self-similar model (©. 
A straight line is drawn between the data points for each cluster for a better visualization of the profiles. The thick black solid line 
represents the median emission measure profile. Right: average gas density profile for the whole sample. The points with errorbars 
represent the median values, while the red points represent the weighted average profile. Note that the errorbars refer to the errors on 
the median emission measure profiles, while the gray shaded region refers to the intrinsic scatter (68.3% confidence level), which amounts 
to ~ 20% 30%) at R 500 (^ 500 )• The long dashed line represents the predictions from the h ydrodynamic numer ical simulations of 

iRoncarelli et all (|20061 ). while the green triangles represents the universal gas density profile from I Eckert et al ] |m3) based on ROSAT 
data. 


Table 2. Gas density, gas density slope /3 = — 1/3 dlog(n e )/dlog(r) and gas fraction at different over densities. The quoted errors refer 
to the errors on the median emission measure profiles. 



A = 500 



A = 200 



A = 100 


sample E 2 n e 

0 

/gas 

E z 2 n e 

0 

/ gas 

E z 2 n e 

0 

/gas 

(10 -4 cm -3 ) 



(10~ 5 cm“ 3 ) 



(10 -3 cm -3 ) 




all 

1.304 

± 

0.034 

0.677 

± 

0.018 

0.130 

± 

0.002 

4.780 


0.199 

0.940 

± 

0.018 

0.150 

± 

0.004 

2.161 

± 

0.130 

1.280 

± 

0.060 

0.158 

± 

0.005 

Ti 

1.310 

± 

0.038 

0.683 

± 

0.017 

0.131 

± 

0.002 

4.744 

± 

0.208 

0.936 

± 

0.018 

0.151 

± 

0.006 

2.116 

± 

0.134 

1.326 

± 

0.060 

0.159 

± 

0.006 

t 2 

1.275 

± 

0.068 

0.646 

± 

0.026 

0.127 

± 

0.003 

5.012 

± 

0.373 

0.945 

± 

0.028 

0.148 

± 

0.008 

1.827 

± 

0.160 

1.375 

± 

0.077 

0.157 

± 

0.008 

t 3 

1.299 

± 

0.095 

0.641 

± 

0.031 

0.128 

± 

0.003 

4.933 

± 

0.473 

0.898 


0.031 

0.150 

± 

0.011 

2.247 

± 

0.186 

1.130 

± 

0.072 

0.161 

± 

0.011 

*1 

1.373 

± 

0.047 

0.642 

± 

0.051 

0.130 

± 

0.002 

4.781 

± 

0.318 

1.004 

± 

0.064 

0.154 

± 

0.005 

2.308 

± 

0.251 

1.121 

± 

0.141 

0.163 

± 

0.008 

Z2 

1.289 

± 

0.074 

0.681 

± 

0.051 

0.133 

± 

0.003 

5.233 

± 

0.344 

1.079 


0.127 

0.156 

± 

0.009 

2.281 

± 

0.214 

1.299 

± 

0.121 

0.161 

± 

0.007 

Z3 

1.287 

± 

0.081 

0.609 

± 

0.029 

0.125 

± 

0.003 

4.639 

± 

0.407 

0.943 

± 

0.036 

0.144 

± 

0.009 

2.144 

± 

0.192 

1.238 

± 

0.069 

0.155 

± 

0.009 

W\ 

1.260 

± 

0.051 

0.657 

± 

0.046 

0.129 


0.002 

4.937 

± 

0.281 

0.922 

± 

0.081 

0.148 

± 

0.006 

2.213 

± 

0.206 

1.212 

± 

0.174 

0.156 

± 

0.006 

W2 

1.324 

± 

0.050 

0.676 

± 

0.019 

0.130 


0.002 

4.805 

± 

0.260 

0.941 

± 

0.018 

0.152 

± 

0.006 

2.261 

± 

0.145 

1.273 

± 

0.058 

0.161 

± 

0.007 

Si 

1.215 

± 

0.038 

0.692 

± 

0.022 

0.125 

± 

0.002 

4.289 

± 

0.246 

1.056 

± 

0.029 

0.141 

± 

0.005 

1.958 

± 

0.162 

1.369 

± 

0.066 

0.148 

± 

0.006 

Si 

1.401 

± 

0.046 

0.625 

± 

0.046 

0.135 

± 

0.002 

5.051 

± 

0.251 

0.853 


0.033 

0.162 

± 

0.004 

2.690 

± 

0.133 

0.985 

± 

0.077 

0.168 

± 

0.006 


Notes on the sample symbols: all : whole sample; T\ : kT < 6 keV; T 2 : 6 < kT < 8 keV; kT > 8 keV; z\ : z < 0.3; Z 2 : 0-3 < z < 0.6; Z 3 : z > 0.6; 
w\ : w < 10 —2 ; W 2 • w > 10 —2 ; f\ : excluding the filament; /2 : including only the filament. 

While the quoted errors refer to the errors on the median, we also report the value of the intrinsic scatter for the whole sample at i?5oo(^200t ^100) : f° r n e 
we have ~ 20%(30%, 38%); for /3 we have ~ 8%(12%, 15%); for /gas we have ~ 15%(25%, 28%). 


© 0000 RAS, MNRAS 000, 000-000 




















12 Morandi et al. 


6.1 Gas density slope and normalization 

Our analysis on these clusters (Figure [Sj reveals that the 
self-similar model is a tenable across a wide redshift and 
mass range. Different from some recent Suzaku results, and 
confirming previous evidence from ROSAT an d Chandra 
llVikhlinin et, al.l ll 9991 . l200fl lEckert et al.1 l2012h , our aver¬ 
age emission measure profile translates into a steep den¬ 
sity profiles. We observe a general trend of steepening in 
the radial profile of the gas density beyond O. 3 R 200 , with 
P ~ 0.940 ± 0.018 at /? 2 oo- This ste epening is in agreemen t 
with previous works by ROSAT (IVikhlinin et al.l 1 199911 . 
but at odds with some recent Suzaku X-ray observations, 
where the electron density decreases steadily with radius, 
approximately following a power-law model wit h /3 ~ 0.7 
llSimionescu et aLll201ll ; iKawaharada et aLll2010ll . 

By using the centroid shift as a proxy of the dynam¬ 
ical state of clusters and the presence of a cool-core (Fig¬ 
ure GO), we did not detect any systematic difference between 
cool-core/relaxed (w < 10~ 2 ) and non-cool core/unrelaxed 
(w > 10 -2 ) clusters at R 500 (Figure [6] on the right). We 
found that NCC clusters do not have on average a higher 
density than CC systems, which might arise due to larger 
clumping factor and /or mer g ing ev ents in disturbed objects. 
In this respect lEckert et al.l (l2012f l noted a clear distinction 
between the two classes in cluster outskirts across all radii, 
out to R 200 , in moderate di sagreeme nt with the present anal¬ 
ysis. As pointed out bv lEckert et al.1 (|2012h . this difference 
in their analysis could be explained by an inaccurate de¬ 
termination of R 200 for unrelaxed clusters in the sample. In¬ 
deed, they use the R 200 —T scaling relation computed under 
the assumption of hydrostatic equilibrium. On the contrary, 
in order to compute R 200 , in the present analysis we used 
the Ey.soo — A /500 relation, which should be nearly insen- 
sitive to differences b etween relaxed and unrelaxed objects 
llKravtsov et, al JfeoQfil i. 

We also do not detect evolution of the normalization 
of the gas density; however, we find hints of evolution of 
the slope of the recovered density profile slope at R 200 is 
P = 1.004 ± 0.064 (1.079 ± 0.127, 0.943 ± 0.036) for 2 < 0.3 
(0.3 < 2 < 0.6, 2 > 0.6, Figure [5] on the left). This differ¬ 
ence can be the result of the evolution of energy feedback 
processes from e.g. AGNs, supernovae and star formation 
as a function of the redshift. Note that current Suzaku and 
ROSAT constraints are limited to very low re dshift, due to 
their poor spatial resolution. In this respect, lEckert et al.1 
(1201214 performed a stacking of the density profiles of a 
sample of ~ 30 clusters at low redshift observed through 
ROSAT, finding P = 0.890 ± 0.026 at R 200 , in agreement 
within 1.7a with our results on the low-redshift subsample. 

A posteriori, our findings on different subsamples sug¬ 
gest that our assumption of self-similarity holds for the ICM 
in the virialization region in clusters, regardless their dynam¬ 
ical states, temperature and redshift. 

Finally, at R 500 the intrinsic scatter of the density 
profiles is ~ 20%, in agreement with previous studies 
JVikhlinin et al.ll2006l : lEckert et al.ll2012l ~). At R 200 , the elec¬ 
tronic density inferred on the whole sample is E~ 2 n e = 
(4.780 ± 0.199) x 10~ 5 cm -3 , with an intrinsic scatter of 
~30%: this value is in agreement with previous ROSAT mea¬ 
sureme nts, i.e. E7 2 n e — (4 .6±0.5) x 10 -5 cm -3 , with a 25% 
scatter dEckert et al.l[2012 h 


6.2 Comparison with simulations 

Comparing our Chandra density profile with numerical sim¬ 
ulations i ncluding cooling, star formation and supernovae 
feedback dRoncarelli et al.]|2006i l. we found that the latter 
fail to reproduce the observed shape of the density pro¬ 
file, predicting density profiles that are significantly too 
steep and with lower normalization compared to the data 
(Figure [5]). This might indicate that non-gravitational ef¬ 
fects are important well outside the core region, whereas 
in the above simulation Roncarelli and collaborators did 
not include AGN feedback and/or preheating. Indeed, re¬ 
cent works dMathews fe Guol 1201 lh indicate that feedback 
mechanisms may be responsible for the deficit of baryons 
in cluster cores, smoothing the accretion pattern and lead¬ 
ing to a flatte r gas distribution. Th is is in agreement wit h 
the findings of lMorandi fc Cul (120141 ') ; lEckert et al.l (l2013al 'l , 
who also argued for an entropy excess significantly beyond 
-R 500 , although we caution the reader that we do not have 
spectroscopic radial temperature/entropy profiles to inde¬ 
pendently validate this scenario. 

The normalization of simulations is also too small: it 
is indeed well known that t hese simulations ove rpredict the 
stellar fraction in clusters dKravtsov fc Borganill2012l ). This 
large reservoir of cold baryons in simulations might po¬ 
tentially explain the discrepancy in the X-ray gas density 
normalization between the present observations and simula¬ 
tions. Consequently, a detailed treatment of gas cooling, star 
formation, AGN feedback, and consideration of gas clump¬ 
ing is required to have realistic models of the outer regions 
of clusters. 

Alternatively, it is possible that X-ray are strongly bi¬ 
ased by the presence of a cold phase and/or clumpy gas 
distribution, or undetected substructures, enhancing the X- 
ray emission and hence the gas density. The CDM scenario 
indeed predicts a picture where clumps at subvirial temper¬ 
ature are infalling along filaments and accreting onto the 
cluster outskirts. Given the low signal-to-noise and/or poor 
spatial resolution of the X-ray/SZ telescope, it might be im¬ 
possible to mask out these substructures e.g. in the X-ray 
analysis, and they might remain undetected. While Roncar¬ 
elli and collaborators eliminated in the simulations the dens¬ 
est clumps (the one per cent of the volume that corresponds 
to the densest SPH particles) that dominate the density and 
surface brightness in the outskirts, this empirical method 
might not fully capture the procedure of masking bright iso¬ 
lated regions from the analysis of observed clusters. This 
procedure indeed depends on the particulars of the observa¬ 
tions, e.g. integration time, X-ray luminosity of the cluster 
and clumps versus the CXB, redshift of the observation, the 
value of the Galactic neutral hydrogen absorption, spatial 
resolution of the X-ray telescope (see i |5.3l and Figure [3] for 
further discussion). 

We then compared our observational results with the 
properties of simulated galaxy clusters produced with th e 
adaptive mesh refinement code ENZ0 dVazza et al.1 120131 / . 
We report that the gas density profiles of simulations is 
roughly in agreement with observations for R > R 500 , the 
former p redicting a slope of ~ 0.97 at the virial radius (see 
lEckert et alj|2012 ). A reasonable agreement with theoretical 
expectations is found when we compare our findings with the 
hydrodynamical numerical simulations of galaxy cluster for- 
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Figure 6. Slope of the gas density 0 = —1/3 dlog(n e )/<flog(r) (left panel) and electron density n e as a function of the radius for different 
populations of the cluster sample. The solid/dashed lines represent the median values for the different subsamples, while the l-cr errors 
for the total sample are represented by the gray shaded region. The panels at the bottom represent the fractional variation of the profiles 
of the subsamples with respect to the whole sample. In order to improve the readability of the figures, we interpolated the data via a 
least-squares constrained spline approximation so as to have continuous functions. 


mation that include radiative cooling, star formatio n , meta l 
enrichment and stellar feedback from iNagai et alj d2007al ). 
These simulations indicate gas density profiles with /3 ~ 0.9 
at f? 2 oo- 

6.3 Unveiling the ongoing large-scale structure 
formation scenario 

Since baryons in the outskirts bear the signature of the con¬ 
tinuous three-dimensional accretion from surrounding fila¬ 
ments, measuring the magnitude of non-equilibrium pro¬ 
cesses such as clumpy gas distribution, the presence of sub¬ 
structures, asphericity in the gas distribution and complex 
accretion patterns provides a neat way to assess the virial- 
ization degree achieved by the ICM. 

In particular, numerical simulations indicate that the 
infall of material into the most massive dark matter haloes 
is not spherical but is expected to be preferentially tun¬ 
neled through the cosmic filaments where the clusters are 
embedded. The cluster mass haloes would indeed acquire 
most of their mass from major mergers along the filaments, 
hence leading to an alignment between the major axis of 
the h ost halo and the large-scale filament dBrunino et al.1 
I2007I) . Therefore an ellipsoidal gas and DM distribution ap¬ 
pears to be a direct outgrow th the ongoing accretion sce¬ 
nario dLimousin et aid 120131 ') . retaining the fingerprints of 
the topology of the cosmic structures. Independent X-ray 
observations confirmed that there is a substantial alignment 
between the major axis of the brightness dist ribution and 
the large-scale filament dMorandi fc Cuil 120141 ). Motivated 


by this cosmic structure formation scenario, we use the po¬ 
sition angle via the moments of the surface brightness dis¬ 
tribution f H4.4I) as a proxy of the direction of the large-scale 
filament where the cluster is embedded. 

In order to gauge the bias implicit in our spherical 
modeling, we repeat the X-ray stacking analysis extract¬ 
ing surface brightness profiles for the 0.7 — 2 keV band for 
two sectors excluding and encompassing the cosmic filament 
(/i and / 2 , respectively). With respect to the azimuthally- 
averaged density profile, the sector /i [fn) shows a den¬ 
sity lower (higher), the net effect being to systematically 
shifting downwards (upwards) the gas density profile, with 
a slope which is flatter when the filament is included (Figure 
©. The difference in gas density normalization between the 
two sectors is ~ 30% at i?ioo- Thus, beyond 7?soo, galaxy 
clusters deviate significantly from spherical symmetry, with 
only small differences between relaxed and disturbed sys¬ 
tems, with the gas located mostly along preferential direc¬ 
tions (i.e., filaments). 

In light of our hypothesis that the elongation of the 
ICM tracks the large-scale filament, the gas density along 
the directions of filaments where the cluster accretes clumpy 
and diffuse materials would then be enhanced in qualitative 
agreement with our analysis. The bias on the X-ray emissiv- 
ity due to subclumps falling into the DM potential well un¬ 
der the pull of gravity should be also more pronounced along 
the direction of the cosmic filaments, if these substructures 
remain undetected in the X-ray anal ysis. In t his res pect, 
as previously pointed out also by I Eckert et al.l (120121) , we 
argued that the shallow density profiles observed in some 
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clusters by Suzaku might be induced by observations in pref¬ 
erential directions (e.g. along filaments), since some these 
observations have been carried along narrow arms, and do 
not reflect the typical behavior of cluster outer regions. 

Finally, we point out that that in Figure [S] the differ¬ 
ences between the /i and /2 sectors originate also within 
J? 500 , in regions where you should naively expect a negligi¬ 
ble impact due to the cosmic filaments. Remembering that 
we recovered the emission measure profiles in circular an¬ 
nuli, this trend in these inner volumes is probably due to 
(mostly) triaxial gas distribution and not to the real pres¬ 
ence of filaments. It is indeed well known that the ICM 
distribution follows isopotential surfaces approximated by 
concentric ellipsoids of decreasing axial ratio towards the 
outer volumes, the dynamics of the gas being driven by the 
gravitational potential well of the DM ( Morandi et al.ll201Cl : 
IMorandi fe Limousinl [201^ ; iMorandi et al.l 20121 ). However, 
in the virialization region (R )> R 500 ) the isopotential sur¬ 
faces of the ICM are expected be nearly spherical, and any 
differences between the /1 and /2 sectors must genuinely 
originate from the impact of the cosmic filaments rather 
from triaxiality. 



0.4 0.6 0.8 1.0 1.2 1.4 

R /^200 


7 GAS FRACTION 

We studied the hot gas mass fraction / gaa (< R) = M gaa (< 
R)/M to t(< R) via stacking of our archival Chandra data. 
We present X-ray measurements of the hot gas fraction out 
to .Rioo by addressing the impact of dynamical state, tem¬ 
perature cluster, evolution with z and asphericity. 

In Figure[7]we present gas fraction profiles for the whole 
sample and for different subsamples. We then compare our 
measure gas fraction with the existing literature out to -R 500 • 
The value of / gaa ,500 = 0.130 ± 0.002 at R 500 . We compare 
our gas fraction to previous studies of / gaa using hydrostatic 
mass estimates. W e used the resu l ts fro m th e samples of 
relaxed clusters bv IVikhlinin et all (120061 1 and lEttori et al .1 
d2009al l. The former measured / gas ,500 = 0.110 ±0.003 for 10 
relaxed clusters observed with Chandra spanning a redshift 
range 2 : = 0.02 — 0.23; the latter / gaa ,500 = 0.106 ± 0.044 via 
52 X-ray luminous galaxy clusters observed with Chandra in 
the redshift range 0.3 — 1.27. These results are in marginal 
disagreement with our determination of the gas fraction at 
-R 500 ■ This stems mostly from the more recent Chandra cali¬ 
brations adopted in the present work, which lead to a lower 
spectral temperature profile (by ~ 10 — 15%) and hence to 
higher values of the gas fraction. Alternatively, our sam¬ 
ple selection criteria could have some ramifications for the 
measurement of / gaa , since we are focusing on luminous 
and massive objects. Indeed the gas fractio n is expec ted to 
be slightly larger for luminous objects (ILandrv et al.ll2013l) . 
while the aforementioned samples includes less massive sys¬ 
tems as well. Our results are also in agreement w ith the 
determination of the gas fraction from I Eckert et all (l2013bl l 
via ROSAT X-ray and Planck Sunyaev-Zeldovich (SZ) ob¬ 
servations (/ gaa ,500 = 0.132 ± 0.005). 

With galaxy clusters at the crossroads of cosmology and 
astrophysics, here we discuss the implications of our findings 
for the use of the gas mass fraction as both cosmological and 
astrophysical probe. 


Figure 7. Cumulative gas fraction / gaa (< R) as a function of 
the radius for different populations of the cluster sample. The 
solid/dashed lines represent the median values for the differ¬ 
ent subsamples, while the 1-cr errors for the total sample are 
represented by the gray shaded region. The two horizzal dot¬ 
ted regions represent, from the top (in dark green) to the bot¬ 
tom (in grey), the f-cr confidence level on the primordial baryon 
fraction from WMAP (12;,/Sl m = 0.167 ± 0.006) and Planck 
(f2j,/n m = 0.155 ± 0.004), respectively. The panel at the bottom 
represents the fractional variation of the profiles of the subsamples 
with respect to the whole sample. In order to improve the read¬ 
ability of the figures, we interpolated the data via a least-squares 
constrained spline approximation so as to have continuous func¬ 
tions. Note that / gaa (< R) scales with the Hubble constant as 
/ gaa oc h~ 3 / 2 . 


7.1 Implications for cosmological studies 

A number of studies have resolved the hot gas mass frac¬ 
tion from X-ray observations in galaxy clusters to place 
constraints on the cosmological parameters, in particular 
Sim, SI a, and the rat io between the pressure and density o f 
the dark energy, w dAllen et al.l l200Sl : lEttori et al.l l2009al ) . 
This hot gas accounts for most of the baryons in clusters: 
the remaining baryons are in stars and intracluster light 
(ICL), and acco unt for a few percent of the total mass 
(ILin et al.l 120031 1. The underlying idea of this cosmologi¬ 
cal test is that massive galaxy clusters are relatively well- 
isolated and gravity-dominated structures, where baryons 
and dark matter have accreted from very large regions of 
~10 comoving Mpc. 

Therefore, their relative baryon budget / gaa should 
be representative of the cosmic value 12(,/fl m = 0.167 ± 
0.006 predicted by the cosmic nucleosynthesis (light ele¬ 
ment formation during the Big Bang) and standard infla- 
tiona ry cosmology, a s inferred from the WMAP 9-year data 
(IBennett et al.l 120131 1. Hence the baryon fraction can be a 
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precise proxy of both the mean dark matter and dark en¬ 
ergy density of the Universe. 

Gas mass fraction measurements are complementary to 
(and competitive with) cosmological constraints using Type 
la supernovae (SNe l a) a nd anisotropi es of the cosmic mi¬ 
crowave background (IA llen et al. 120111) . However, they rely 
on the key assumptions that the hot gas fraction within a 
given radius is independent of redshift; can be converted to 
the universal fraction by correction factors which measure 
the “baryon depletion “ from cosmological simulations and 
the fraction of baryons in stars and galaxies (not directly 
constrained via X-ray observations); and can be measured 
reliably. 

In order to understand how biases on the hot gas frac¬ 
tion / g aa, stellar fraction /stars and depletion parameter b 
(the parameter used to convert the hot gas fraction to the 
cosmic ratio) propagate into the determination of, e.g. fl m , 
we can reverse-engineer the total matter content using the 
observed baryon fraction in clusters: 

i 1 m — b Hb/(/gas T /stars ). (10) 

with /stars ~ 0.12 dGonzalez et al.|[2013l ). 

After correcting the hot gas fraction for the baryons in 
stars, our reconstructed average baryon fraction increases 
with radius and reaches the value /b,ioo = 0.177 ± 0.005, 
slightly larger than the primordial value. 

We point out that in our stacking analysis we assume a 
priori that the gas fraction, and hence the gas densities, do 
not depend on the system temperature (Equation [3]). This 
hypothesis is verified a posteriori in light of a lack of de¬ 
pendency of the gas fraction on the cluster temperatures. 
This is probably explained by the relatively narrow temper¬ 
ature range spanned in our sample (> 3 keV), such that the 
clusters in our sample should show little dependence of the 
physical properties on the global temperature. Indeed, while 
the gas fraction in groups and intermediate mass clusters im¬ 
plies an increasing trend of th e cu mulative gas mass fraction 
with temperature (iSun et al.ll2009l ). it is less clear whether 
such a trend persists for rel a tively mas sive objects . Mea - 
surements bv IVikhlinin et all (120061 ) and I Allen et al] (120081 ) 
are both consistent with being constant with cluster tem¬ 
perature, for objects with kT > 5 keV, in agreement with 
our findings. Note that most of the clusters (264 out of 320) 
in our sample have temperatures greater than 5 keV. This 
trend is also confirmed by the level of the intrinsic scatter 
of /gas (~ 15% at J? 5 oo, ~ 25% at R 200 ), which does not 
appear to depend on the selected subsamples. 

Our results also demonstrates that, when computed 
within the whole cluster virialization region, there is no rel¬ 
evant dependency of the gas fraction on cluster dynamical 
state at radii larger than R 500 . These re sults are in tension 
with the findings of lEckert et al.l (l2013bl) . who finds signifi¬ 
cant differences between the baryon fraction of relaxed, cool- 
core (CC) systems and unrelaxed, non-cool core (NCC) clus¬ 
ters out to R 200 . This difference might arise since they ex¬ 
plicitly use hydrostatic equilibrium (by combining the pres¬ 
sure from SZ Planck data with the gas density from ROSAT) 
to calculate total masses, which can bias upward the gas 
mass especially in unr elaxed systems llBattaglia et al.ll2013l ; 
iRoncarelli et al.ll2013f) . 

Appreciably high / gas (larger than the cosmic value) 
towards the virial radius have been also reported in re¬ 


cent studies based on Suzaku data (e.g., ISimionescu et all 
H), calling for a significant clumpy distribution of the 
ICM and the presence of non-thermal pressure support, in 
moderate disagreement with the present analysis. In this 
respect, understanding the unbiased value of the baryonic 
budget in clusters has clearly important ramifications in 
the use the hot gas fraction as a potential stand a rd can - 
dle for the cosmology. For example, lEttori et alj J2009ah 
found a tension between the inferred value of (jointly 
constrained with Ha) via the hot gas fraction in X-ray clus¬ 
ters (fi m = 0.32 /q'q 2) and the value from WMAP 5-year 
data (0.279 ±0.013), with the former exceeding the latter. If 
the WMAP best-fit results are assumed to fix the cosmolog¬ 
ical parameters, this tension is expected to arise from biases 
due to non-thermal pressure support and ICM dumpiness. 

We finally observe no evolution of the gas fraction with 
redshift. Measurements of the apparent lack of evolution of 
the cluster X-ray gas mass fraction has been used to probe 
the acceleration of the Universe, through the d ependence of 
/gas on the assumed distances to the clusters dEttori et al.l 
120061 . l2009al ) . This assumption has never been tested obser- 
vationally. Our results suggest, for the first time, that this 
assumption in tenable for cosmological purposes. 

7.2 Implications for the cluster astrophysics 

The nearly closed-box nature of deep potentials of massive 
galaxy clusters makes them ideal laboratories to study non- 
gravitational processes, in particular galaxy feedback from 
supernovae and AGNs, operating during galaxy formation 
and their effects on the surrounding intergalactic medium. In 
the literature there has been a long-standing debate on the 
baryon content of clusters, where most of the studies indicate 
that the cluster baryon budget is typically lower than the 
cosm i c baryon fraction llVikhlinin et al.l 120061 : I Allen et al.l 
l200Sl : lEttori et al.ll2009al ). although other analyses find that 
the b aryon content is consistent with the cosmic ratio at 
Rsoo dMiller et al .11201 ‘J : Ibandrv et al.ll2013l ). A possible ex¬ 
planation of the observed deficit of baryons was that non- 
gravitational processes associated with galaxy formation, 
feedback from supernovae, AGN, star formation, or galac¬ 
tic winds may push t he ICM towards the clu s ter outskirts 
(R ~ R200, see, e.g. Metzler fe Evrardlll994l ; iBialek et aid 
l200ll ; iMcCarthv et al. 20071 : Sun et ak[ 20091 ). Nevertheless, 
we stress that current studies are mostly limited to the in¬ 
ner volumes of the clusters (R < R2500 — R500), because of 
the very low surface brightness of the X-ray signal in the 
outskirts. 

As previously pointed out, we do not observe evolution 
of the gas fraction with redshift and an overall baryon bud¬ 
get slightly larger than the cosmological value within i?ioo- 
This might dispel questions which arose in the literature on 
the whereabouts of these putative missing baryons. This also 
confirms that the most massive clusters, which are relatively 
well-isolated and gravity-dominated structures, should re¬ 
tain all their gaseous matter. Moreover, while any feedback 
process will be more significant in the inner regions of low- 
mas s systems when compared to the binding energy of the 
gas dSun et al.[|2009l ). we checked that in the outskirts there 
seems to be no particular dependency on the cluster temper¬ 
ature or dynamical state. A posteriori, these findings suggest 
that our assumption of self-similarity holds for the ICM in 
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the virialization region in clusters, regardless their dynami¬ 
cal states and mass. 

We point out that an additional baryonic component 
might not be detected, raising the possibility that a non- 
negligible fraction of baryonic mass may be hidden in a 
form that is difficult to measure, such as cool gas. As we 
observed in 1 16.21 the presence of dense and cold clumps at 
temperatures kT <, 1 keV in the ICM might remain un¬ 
detected, since their exponential bremsstralilung cutoff will 
be at energies ~ kT/( 1 + z) below the energy threshold of 
the X-ray detector (~ 0.5 keV for Chandra). Moreover, if 
the multiphase structure of the ICM turns out to be im por¬ 
tant as numerical simulations suggest (see, e.g.. iRasia et al.l 
l2014h . we might actually underestimate baryons due to the 
single-temp erature m odelling of the ICM 1 46.211 . In this re¬ 
spect, lAfehordFet^I] (j 2007t f stacked the SZ effect signals of 
WMAP observations of a large sample of massive clusters to 
constrain the thermal energy of the ICM. With the aid of 
hydrodynamic simulations, they converted this thermal en¬ 
ergy into an estimate of the hot gas fraction of clusters, and 
found /gas ,200 = 0.109±0.013. This implies a deficit (~ 35%) 
of baryons from the hot ICM (without correcting for /stars), 
calling for a cool phase of the ICM. This effect could be 
competitive with the enhancement of the gas fraction due to 
inh omogeneities of the gas distribution (~ C 0 ' 5 ~ 1.2 — 1.4, 
see iRoncarelli et al.l ll2013ll ; iBattaglia et al.l J2013I ') and be¬ 
low discussion). 

We also observe that the gas faction along the direc¬ 
tion of the filaments is systematically higher (~ 10 — 15%) 
than the remaining volume. This picture is supported by our 
analysis of azimuthal variations of the gas density in cluster 
outskirts (» which suggests that even CC clusters ex¬ 
hibit significant departures from spherical symmetry around 
i? 200 - Consequently, a full azimuthal coverage is required to 
study the baryon budget of cluster outer regions. 

Next, we investigated the biases on our gas fraction 
measurements due to the effects of inhomogeneities of the 
gas distribution (’’dumpiness”) and non-thermal pressure. 
For the former bias, observations and simulations show that 
the gas clumping factor C lies in the range ~ 1.3 — 2 at R 200 
llNagai fe Laul 1 20 111 : I Morandi et al.l l2013bl : iMorandi fe Cuil 
120141 ), which wil l bias upwards the gas mass fraction of about 
C 0,5 ~ 1.2—1.4 (IBattaglia et al]|2013l;l Roncarelli et al.ll2013l : 


IMorandi et al.l l2013bl : Morandi fe Cuil 2014 ). For the latter 
bias, our measurements are also affected by violation of hy¬ 
drostatic equilibrium 1 45.211 . A bias b ~ 0.1 translates into 
a bias (upwards) of ~ 14% of the gas fraction at Raoo- We 
stress that the real challenge in future gas fraction measure¬ 
ments will be to account for these systematics, since nowa¬ 
days X-ray studies customarily hinge on the assumptions of 
hydrostatic equilibrium, single-temperature spectral models 
and homogeneous distribution of the gas density. If these 
estimates of the biases due to non-thermal pressure and in 
particular dumpiness are correct, we might non-negligibly 
overestimate the gas fraction, leading to an actual deficit of 
baryons in clusters with respect to the cosmic value. How¬ 
ever, this deficit might be dispelled if we change the cosmol¬ 
ogy from WMAP to Planck, as we discuss below. 

Finally, we observe that the differences between the 
WMAP 9-year and our cosmological parameters affect 
the estimates of the gas mass fraction through the vari¬ 
ation of the angular diameter distance d ang (/gas oc 


d a , U g(Ho,fl ln ,flA,w) 3 / 2 ). If we assume the WMAP 9-year 
parameters, high -z clusters will have a slightly larger (<^ 2 
percent) gas fraction. 

Moreover, we analyzed the effects of changing the cos¬ 
mology from WMAP to Planck. Plank data seem to sug¬ 
gest an higher value of the total matter density of the 
Universe (> 0.3) and a lower value of the Hu bb le con¬ 
stant with respect to WMAP dPlanck Collaboration et al.l 
I2013T) . This translates into a lower value of the cosmic 
baryon fraction flb/fim = 0.155 ± 0.003. On the other 
hand, the net effect of switching from WMAP to Planck 
cosmology is to increase the measured gas fractions by 
(dang,Plank/dang, wmap ) 3/2 ~ 1-05 for high-z systems, and 
stellar fractions by (da ng ,piank/dang,WMAp) ~ 1-03. Clearly, 
the most striking effect of the use of the Planck cosmologi¬ 
cal para meters is to increase t he baryon content within Rioo 
(see also lGonzalez et al]|2013h . 


8 SUMMARY AND CONCLUSIONS 

Cluster outskirts are the regions where the transition be¬ 
tween virialized cluster gas and accreting material from the 
large-scale structure occurs. In this paper, we have presented 
our analysis of a sample of 320 clusters (z = 0.056 — 1.24) 
from the Chandra archive, focusing on the properties of the 
gas in cluster outskirts. 

We have exploited the large archival dataset of Chan¬ 
dra clusters to trace the gas density and gas fraction of the 
intracluster gas out to Rioo- We measured the median of the 
distribution of the (renormalized) emission measure profiles 
EM(R/R 200 ) to detect signal in the outer volumes and mea¬ 
sured the typical gas density, gas slope and gas fraction. We 
investigated the evolution of the physical properties as a 
function of dynamical state, redshift, cluster temperatures, 
as well as their azimuthal variation. Finally, we compared 
our average density and gas fraction with the results of nu¬ 
merical simulations. 

In tension with some recent Suzaku results, and confirm¬ 
ing previous evidence from ROSAT and Chandra, we observe 
a steepening of the density profiles beyond Rsoo- We report 
that galaxy clusters deviate signihcantly from spherical sym¬ 
metry, but with negligible differences between relaxed and 
disturbed systems. Thus, our findings on different subsam¬ 
ples suggest that our assumption of self-similarity holds for 
the ICM in the virialization region in clusters, regardless 
their dynamical states, temperature and redshift. 

Comparing our results with numerical simulations, we 
find that the latter simulations fail to reproduce the gas 
distribution, even well outside cluster cores. Simulations in¬ 
cluding cooling and star formation convert a large amount 
of gas into stars, which results in a low gas density and gas 
mass fraction with respect to the observations. 

We then measured the hot gas fraction in galaxy clus¬ 
ters observed with Chandra out to Rioo- After converting 
the hot gas fraction to the total baryon budget in clusters, 
we measured a value of the baryon fraction which slightly ex¬ 
ceeds the cosmic baryon fraction. However, a careful analysis 
of the systematics due to non-thermal pressure and dumpi¬ 
ness suggest that we might non-negligibly overestimate the 
total baryon budget, leading to an actual deheit of baryons 
in clusters with respect to the cosmic value. 
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This study has important implications on our knowl¬ 
edge of the state of the gas in cluster outskirts, and promises 
to help advance cosmological studies based on the cluster 
growth over cosmic time in the coming decade. 


9 FUTURE WORKS 

Our work opened the way to a novel method, the stacking 
the X-ray signal of cluster outskirts, exploiting the high level 
of self-similarity of the cluster outer volumes. The total ex¬ 
posure time is ~ 20 Ms for the whole sample, which is one 
order of magnitude larger than the deepest Chandra obser¬ 
vations (e.g. compared to the 2 Ms of A133, see Vikhlinin 
et al. in prep.), and 2-3 order of magnitudes larger than an 
average X-ray observation. Clearly, the most striking rami¬ 
fication of this analysis is the ability to greatly exploit the 
Chandra database, while no single observation could achieve 
the same scientific goals. 

In future papers will use the stacking approach to infer 
constraints on the cosmological parameters via the analysis 
of the X-ray in the outskirts (Morandi et al. in prep.). This 
novel approach will pin down the cosmological parameters, 
in particular fl m and w (the ratio between the pressure and 
density of the dark energy), with a greatly improved accu¬ 
racy with respect to the current analyses of clusters, char¬ 
acterized by relatively smaller samples and limited to the 
interiors of the sources (e.g. out to -R 2500 , As 00 )• We plan to 
use the X-ray stacked data to improve further our under¬ 
standing of how the integrated Compton parameter scales 
with total cluster mass and the corresponding Compton SZ 
parameter from Planck , when available; and to constrain the 
temperature and the metallicity profiles in the outer regions, 
in order to understand the thermodynamic of the outer re¬ 
gions. 
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APPENDIX A: DISTRIBUTION OF THE 
X-RAY COUNTS IN THE SAMPLE 

We present the number of net X-ray counts per annulus for 
our sample. We also show the distribution of the number of 
net X-ray counts for whole dataset. 


APPENDIX B: SPECTRAL TEMPERATURES 
DEPENDENCY ON ATOMIC DATABASE AND 
ABSORPTION 

We present a comparison between spectral temperatures 
and metallicities recovered via the AtomDB databases 2.0.2 
and 1.3.1, and via different models of X-ray absorption by 
the ISM (Figure IBII) . Note that in our work w e adopted 
the APEC emissivity model dFoster et alj 120121 1 and the 
AtomDB (version 2.0.2) database of atomic data, and we 
emplo yed the solar abundance ratios from lAsplund et al.l 
(120091b We also used the Tuebingen-Boulder absorption 
model (tbabs) for X-ray absorption by the ISM. 

Finally, we discuss the bias due to our choice to fix 
the hydrogen column density Nh to the Galactic value 
by using the Leide n/Argentine/Bonn (LAB) Hl-survey 
dKalberla et al.l 120051 1. This radio surveys measure the hy¬ 
drogen column density from the HI-21 cm line, only pro¬ 
viding the neutral hydrogen along the line of sight, while 
the molecular a n d ion ized hydrogen is not accounted for. 
IWillingale et alj (|2013l ) provided a method to account for 
the molecular hydrogen via the dust extinction measured in 
the B and V bands, and by using X-ray afterglows of gamma 
ray bursts. Although the neutral hydrogen usually con¬ 
tributes most of the total hydrogen for Nh <, 10 21 cm~ 2 , the 
impact of molecular hydrogen should be important for larger 
column densities. For each cluster, we performed mock simu¬ 
lations of absorbed spectra by assuming the our average gas 
density p rofile ( jQIl ), the average temperature profile as mea¬ 
sured bv lvikhlinin et abl d200al and hydrogen column den¬ 
sity values including both neutral (from the aforementioned 
radio surveys) and m olecular hydrogen (as determined by 
IWillingale et al.ll2013l b These mock spectra have been then 
fit by fixing the hydrogen column d ensity to the G alacti c 
value as determined by the Hl-survey dKalberla et al.|[2005l 4 . 
For Nh < 6 X 10 20 cm -2 , the bias (upwards) on the recov¬ 
ered spectral temperature is <) 1 — 2%, which translates into 
a bias of similar magnitude (downwards) on the gas density 
n e , 200 - However, for Nh > 10 21 cm -2 the bias is apprecia¬ 
ble (j> 10%). Note that most of the clusters (290 out of 320) 
have column densities smaller than 6 x 10 2O cm -2 , while only 
3 clusters have Nh > 10 21 cm -2 . Hence the overall system- 
atics on the gas density recovered via stacking is estimated 
< 2% at f? 200 - 


APPENDIX C: DEFINITION OF COOLING 
CORE SYSTEMS 

We classify clusters as cool cores on the basis of density 
Ep 2 n e fi at 0.03/?5oo- The threshold we use to define a cool 
core systems is FT 2 n e ,o > 1.5x 10~ 2 cm~ 3 . We compared ou r 
definition of cool core systems with lCavagnolo et al] (l2009h . 
whose sample shares 93 objects in common with ours. They 
use the entropy threshold A'o = 30 — 50 keV cm 2 to approx¬ 
imately demarcate the division between CC and NCC. In 
Figure ED we present a comparison between the two defi¬ 
nitions of cooling core cluster, along with the definition of 
cool core based on the cooling time. As shown in Figure ED 
our de finition of cool cores is r oughly agrees with the defini¬ 
tion of lCavagnolo et abl (l2009h , while the definition based on 
cooling time falls in between. In the end, our conclusions on 


© 0000 RAS, MNRAS 000, 000-000 






























18 Morandi et al. 


R500 ^100 

0.5 


0.4 


U 

o 

CD 

1 0.2 

o 


0.1 


0.0 

0.1 1.0 

^/^200 

Figure Al. Left panel: number of net X-ray counts per annulus in the energy range 0.7-2.0 keV for the 320 clusters of our sample. The 
shaded gray band represents the measured total background in each annulus. Central panel: distribution of total net X-ray counts for 
the 320 clusters used in this work. For all the panels blue and red curves (dots) correspond to clusters classified as cool core and non-cool 
core, respectively (see rot . 





1 1 1 I—f 



l 

nH-0.01(x10 22 ) 

. nH=0.04(x10 22 ) 

- nH=0.10(x10 22 ) 

- nH=1.00(x10 22 ) 

' 

i t 

t 1 

\ 

\ 

\ 

\ L 


\ 

■I 

V\ 

1 

K 1 

1 \l 

1 

' \ 

\ 


_ 

\ 

\ 

f- 1 

1 


\ A :• J' 

\ w \/\ A /- 

\ / v v \S V A ' 

V J\ I 



^ 


1 

T (keV) 


10 


Figure Bl. Left panel: comparison between the spectral temperatures and metallicities recovered via the AtomDB databases 2.0.2 and 
1.3.1 via mock data. Right panel: comparison between the spectral temperatures recovered via different models of X-ray absorption by 
the ISM via mock data (tbabs v.s. phabs). 
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Figure Cl. Comparison between the central entropy Kp from 
the su bsample of 93 objects in common with ICavagnolo et al. 
(120091) with our gas density at 0.03-R500- The threshold we use 
to de fine a cool core s ystem s is E^ 2 n e ,o > 1.5 x 10 -2 cm~ 3 , 
while ICavagnolo et all (120091 ) use the entropy threshold Kq ~ 
40 keV cm 2 to approximately demarcate the division between CC 
(cool core clusters) and NCC (non-cool core clusters). The blue 
line divide the plot in two region with cooling time greater (upper 
region) and smaller (lower region) than 0.4 Hubble times. 


the CC/NCC classification are not sensitive to our definition 
of cool cores. 


APPENDIX D: COMPARISON OF Mgas ,500 
WITH THE LITERATURE 

We compared our gas mass at R 500 (recovered via deprojec¬ 
tion of individual e missio n measure profiles) with the results 
of lMaughan et al] (I2008I) . who se samp l e shares 95 o bjects in 
common with ours, and with iMartino et al] (120141 ). whose 
sample contains 32 objects in common with us (Figure TD II) . 
We found no significant differences between the recovered 
gas masses m, with a scatter of ~ 5%. We performed the 
comparison at the same R 500 from the literature, and the gas 
mass has been recovered from individual gas density profile 
rather than the stacked profiles. 
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The galaxy cluster outskirts probed by Chandra 21 


Table 1. The X-ray properties of the galaxy clusters in our sample. For each object different columns report the cluster name, the 
redshift z, the exposure time t e xp, the identification number of the Chandra observations, the spectroscopic temperature T ew in the 
radial range 0.15 — 0.75 # 500 , the centroid shift (in units of -R 500 ), and a flag for the presence or not of a cooling core (labeled CC and 
NCC, respectively). 


Cluster 

z 

texp 

ks 

ID 

T e „ 

(keV) 

w 

(xl0“ 2 ) 

CC/NCC 

A85 

0.056 

272.4 

904 4881 4882 4883 4884 4885 4886 4887 

4888 15173 16264 15174 16263 

6.45 + 0.86 

0.49 + 0.18 

CC 

A133 

0.057 

2393.0 

13442 13443 13444 13445 13446 13447 13448 13449 

13450 13451 13452 13453 13454 13455 13456 13457 
14333 14338 14343 14345 14346 14347 14354 13391 

13392 13518 3183 3710 12177 12178 12179 9897 

3.76 ± 0.25 

0.26 ± 0.09 

CC 

A644 

0.070 

48.9 

2211 10420 10421 10422 10423 

8.12 + 1.12 

1.56 ± 0.57 

CC 

A401 

0.075 

153.9 

10416 10417 10418 10419 14024 

5.44 + 0.86 

1.03 ± 0.38 

NCC 

A2029 

0.076 

29.2 

6101 10434 10435 10436 10437 

7.38 ± 0.92 

0.21 ± 0.08 

CC 

A1650 

0.084 

218.4 

7691 5822 5823 6356 6357 6358 7242 10424 

10425 10426 10427 

5.89 + 0.86 

0.38 ± 0.14 

CC 

A1068 

0.137 

19.5 

13595 13596 13597 13598 

6.38 + 1.26 

0.29 ± 0.11 

CC 

A2276 

0.141 

39.5 

10411 

3.40 ± 0.48 

0.49 ± 0.18 

NCC 

A1413 

0.143 

155.6 

5003 12194 12195 12196 13128 1661 5002 537 

7696 

7.78 ± 0.22 

0.57 + 0.21 

CC 

A3402 

0.146 

19.8 

12267 

3.09 + 0.55 

1.39 ± 0.51 

NCC 

A2409 

0.148 

10.2 

3247 

5.40 ± 0.40 

1.09 ± 0.40 

NCC 

RXCJ0449.9-4440 

0.150 

19.8 

9417 

4.23 + 0.47 

0.66 ± 0.24 

NCC 

A2204 

0.152 

106.5 

7940 6104 12895 12896 12897 12898 

10.38 + 0.14 

0.18 ± 0.06 

CC 

A907 

0.153 

99.0 

3185 3205 535 

6.60 ± 0.24 

0.17 + 0.06 

CC 

RXJ1000.4+4409 

0.154 

18.5 

9421 

3.42 ± 0.41 

3.42 ± 1.25 

CC 

A3866 

0.154 

9.0 

15162 

6.01 + 0.84 

0.15 ± 0.05 

CC 

Zw8284 

0.156 

22.0 

15118 

5.00 + 0.46 

0.50 + 0.18 

NCC 

RXCJ0528.2-2942 

0.158 

19.3 

9418 

6.41 + 0.48 

0.75 ± 0.27 

NCC 

A1240 

0.159 

51.4 

4961 

6.85 + 0.85 

2.05 ± 0.75 

NCC 

A2259 

0.164 

10.0 

3245 

6.37 + 0.44 

0.60 ± 0.22 

NCC 

A2445 

0.165 

26.7 

12249 

5.62 + 0.77 

0.63 ± 0.23 

NCC 

A853 

0.166 

24.7 

12250 

6.50 + 1.27 

0.52 ± 0.19 

NCC 

A1201 

0.169 

47.4 

9616 

6.47 + 0.41 

2.54 ± 0.93 

NCC 

Z808 

0.169 

18.8 

12253 

6.95 + 1.65 

0.24 ± 0.09 

CC 

ZwC115 

0.170 

26.7 

12251 

6.11 + 0.63 

0.43 ± 0.16 

NCC 

A1204 

0.171 

23.3 

2205 

4.32 + 0.39 

0.22 ± 0.08 

CC 

A1914 

0.171 

38.6 

3593 12197 12892 12893 12894 

9.18 + 0.91 

0.64 ± 0.23 

CC 

A586 

0.171 

9.9 

11723 

6.85 + 1.01 

0.48 ± 0.18 

CC 

RXCJ0616.3-2156 

0.171 

23.8 

15100 

6.73 + 0.35 

2.66 ± 0.97 

NCC 

RXJ1750.2+3505 

0.171 

19.8 

12252 

6.14 + 1.01 

0.40 + 0.15 

CC 

A3140 

0.173 

20.0 

9416 

6.25 + 0.81 

0.47 + 0.17 

NCC 

A2218 

0.176 

46.3 

7698 1666 

7.96 + 0.45 

1.00 + 0.36 

NCC 

MS0906.5+1110 

0.180 

34.5 

924 7699 

5.57 + 0.44 

0.99 ± 0.36 

NCC 

A665 

0.182 

102.3 

12286 13201 15147 15148 

8.05 + 0.23 

4.15 ± 1.51 

NCC 

A2187 

0.183 

18.3 

9422 

7.55 + 0.74 

1.61 ± 0.59 

NCC 

A1689 

0.183 

196.7 

6930 7289 7701 540 1663 5004 

10.27 + 0.26 

0.52 ± 0.19 

CC 

A0598 

0.186 

19.8 

10442 

6.28 + 1.39 

0.78 ± 0.29 

CC 

A383 

0.187 

29.2 

2320 524 

5.47 + 0.32 

0.20 + 0.07 

CC 

RXCJ0331.1-2100 

0.188 

20.0 

10790 9415 

6.03 ± 0.73 

0.80 ± 0.29 

CC 

CIZAJ2242.8+5301 

0.192 

195.9 

14019 14020 

9.37 + 0.39 

2.99 ± 1.09 

NCC 

A4023 

0.193 

22.8 

15124 

7.91 ± 1.12 

2.56 ± 0.94 

NCC 

MS0839.8+2938 

0.194 

4.8 

7702 

4.49 ± 0.50 

0.91 ± 0.33 

CC 

A115 

0.195 

360.2 

13458 13459 15578 15581 3233 

7.85 ± 0.22 

4.12 ± 1.50 

CC 

A2507 

0.196 

33.3 

12248 

9.86 + 1.62 

0.53 ± 0.19 

NCC 

RXJ2247+0337 

0.200 

49.0 

911 

3.04 + 0.47 

1.46 ± 0.53 

NCC 

A3322 

0.200 

15.6 

15111 

7.47 + 1.05 

0.54 ± 0.20 

NCC 

A520 

0.202 

408.1 

9424 4215 9426 9430 528 

8.74 + 0.25 

1.88 + 0.68 

NCC 

A3399 

0.203 

24.5 

15125 

6.76 ± 0.76 

1.41 ± 0.51 

NCC 

ZwC11829.3+6912 

0.204 

64.0 

10412 10931 

4.85 ± 0.64 

0.83 ± 0.30 

NCC 

A963 

0.206 

5.1 

7704 

7.40 + 1.75 

0.38 + 0.14 

CC 

A209 

0.206 

20.0 

522 3579 

8.69 ± 0.64 

0.76 ± 0.28 

NCC 

A223 

0.207 

44.8 

4967 

4.16 + 0.43 

1.28 ± 0.47 

NCC 

RXJ0439+0520 

0.208 

38.1 

527 9369 9761 

5.14 + 0.55 

0.38 ± 0.14 

CC 
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Table 1. Continued. 


Cluster 

z 

^•exp 

ks 

ID 

T ew 

(keV) 

w 

(xlO -2 ) 

CC/NCC 

MACSJ0547.0-39040.210 

21.5 

3273 

4.08 ± 0.39 

1.77 + 0.65 

cc 

ZWCL2701 

0.210 

5.1 

7706 

4.70 + 0.62 

1.07 + 0.39 

cc 

G286.58-31.25 

0.210 

21.9 

15115 

7.03 ± 0.95 

1.42 ± 0.52 

NCC 

A1430 

0.211 

21.8 

15119 

6.74 ± 0.73 

2.06 ± 0.75 

NCC 

A222 

0.213 

45.1 

4967 

3.88 ± 0.30 

3.65 ± 1.33 

NCC 

A1246 

0.213 

5.0 

11770 

6.03 + 0.99 

1.08 ± 0.40 

NCC 

A1423 

0.214 

35.5 

538 11724 

6.61 + 0.53 

0.67 ± 0.24 

CC 

MS0735.6+7421 

0.216 

475.3 

10468 10469 10470 10471 10822 10918 10922 

6.64 + 0.14 

0.47 + 0.17 

CC 

A773 

0.217 

60.2 

5006 3588 533 13591 13592 13593 13594 

7.53 ± 0.49 

1.05 ± 0.38 

NCC 

RXCJ1947.3-7623 

0.217 

19.8 

15102 

7.16 + 1.15 

0.83 ± 0.30 

CC 

A3084 

0.219 

19.8 

9413 

6.46 ± 0.80 

0.78 ± 0.29 

NCC 

A3017 

0.220 

14.9 

15110 

7.11 ± 1.14 

1.61 ± 0.59 

CC 

RXCJ0510.7-0801 

0.220 

20.7 

14011 

6.44 ± 0.43 

1.71 ± 0.63 

NCC 

A368 

0.220 

18.4 

9412 

8.28 + 1.25 

0.49 ± 0.18 

CC 

MS1006.0+1202 

0.221 

67.5 

925 13390 

6.30 ± 0.48 

2.16 ± 0.79 

NCC 

AS0592 

0.222 

19.9 

9420 

8.91 ± 1.08 

0.72 ± 0.26 

CC 

RXCJ1514.9-1523 

0.223 

59.0 

15175 

9.41 + 0.69 

1.83 ± 0.67 

NCC 

A1763 

0.223 

19.6 

3591 

8.63 ± 0.82 

1.16 ± 0.42 

NCC 

A1942 

0.224 

60.6 

3290 7707 

4.92 ± 0.46 

1.14 + 0.42 

NCC 

A2261 

0.224 

24.3 

5007 

8.09 ± 0.62 

0.28 + 0.10 

CC 

A1895 

0.225 

18.7 

15129 

4.93 ± 0.84 

1.26 ± 0.46 

NCC 

A2219 

0.226 

152.5 

13988 14355 14356 14431 14451 7892 

10.81 + 0.31 

1.23 ± 0.45 

NCC 

ZwC10823.2+0425 

0.226 

21.4 

10441 

4.83 + 0.51 

0.57 + 0.21 

NCC 

RXCJ0118.1-2658 

0.228 

19.7 

9429 

7.45 ± 0.57 

2.38 ± 0.87 

NCC 

RXJ0220.9-3829 

0.229 

19.9 

9411 

5.10 + 0.65 

0.37 + 0.13 

CC 

RXJ1234.2+0947 

0.229 

29.3 

11727 539 

6.30 + 1.01 

1.99 ± 0.73 

NCC 

A141 

0.230 

19.9 

9410 

6.01 + 0.81 

1.24 + 0.45 

NCC 

A267 

0.230 

19.9 

3580 

7.92 ± 0.62 

1.46 ± 0.53 

NCC 

A2111 

0.230 

31.2 

544 11726 

8.66 + 0.75 

2.51 ± 0.92 

NCC 

A2667 

0.230 

19.4 

13599 13600 13601 13602 

6.31 ± 1.05 

0.34 + 0.13 

CC 

RXJ0439.0+0715 

0.230 

18.8 

3583 

6.30 + 0.62 

0.90 + 0.33 

CC 

A1682 

0.234 

29.4 

3244 11725 

6.47 + 0.70 

0.85 ± 0.31 

NCC 

A2146 

0.234 

375.5 

12245 13120 12246 13138 12247 13023 13020 13021 

7.28 + 0.20 

3.63 ± 1.32 

CC 

ZwC12089 

0.235 

9.1 

7897 

5.13 + 0.71 

0.23 + 0.08 

CC 

RXJ2129.6+0005 

0.235 

39.6 

552 9370 

7.78 ± 0.40 

0.71 ± 0.26 

CC 

A2465 

0.245 

69.2 

14010 15547 

3.97 + 0.39 

1.53 ± 0.56 

NCC 

A2125 

0.247 

117.3 

2207 6891 7708 

3.56 + 0.25 

2.53 ± 0.92 

NCC 

A2485 

0.248 

19.8 

10439 

5.86 + 0.53 

0.35 ± 0.13 

NCC 

A2645 

0.251 

19.0 

14013 

6.35 + 1.01 

0.83 ± 0.30 

NCC 

A1835 

0.253 

192.4 

6880 6881 7370 

10.00 + 0.29 

0.35 ± 0.13 

CC 

A521 

0.253 

88.4 

12880 13190 

7.51 + 0.41 

5.95 ± 2.17 

NCC 

A3088 

0.253 

18.9 

9414 

8.22 + 1.04 

0.28 ± 0.10 

CC 

A68 

0.255 

10.0 

3250 

8.97+1.91 

1.14 + 0.42 

NCC 

MS1455.0+2232 

0.259 

108.0 

4192 7709 543 

5.26 + 0.22 

0.32 ± 0.12 

CC 

RXCJ1327+0211 

0.262 

6.7 

11771 

8.76 + 1.70 

1.10 + 0.40 

NCC 

RXJ0528.9-3927 

0.263 

109.9 

4994 15658 15177 

10.64 + 0.63 

1.59 ± 0.58 

CC 

Z5768 

0.266 

37.3 

14014 7898 

3.29 ± 0.40 

3.16 ± 1.15 

NCC 

MACSJ2211.7-03490.270 

17.7 

3284 

11.97 + 2.69 

0.86 ± 0.32 

CC 

RXCJ0532.9-3701 

0.271 

24.8 

15112 

9.20 + 1.24 

0.40 + 0.15 

CC 

RXCJ0303.7-7752 

0.274 

33.6 

15113 

9.83 ± 0.74 

1.37 + 0.50 

NCC 

SDSSJ1233+1511 

0.275 

12.1 

11761 

5.17 + 0.61 

3.46 ± 1.26 

NCC 

A1622 

0.275 

12.9 

11763 

4.76 + 0.87 

1.24 + 0.45 

NCC 

RXJ0142.0+2131 

0.277 

19.4 

10440 

8.50 + 0.66 

1.21 ± 0.44 

NCC 

A2631 

0.278 

26.0 

3248 11728 

7.15 + 0.82 

1.84 + 0.67 

NCC 

ACTJ0235-5121 

0.278 

19.8 

12262 

8.18 + 0.76 

2.28 + 0.83 

NCC 

SDSSJ0104+0003 

0.278 

15.0 

11765 

6.49 ± 0.80 

2.83 ± 1.03 

NCC 

A1576 

0.279 

43.4 

7938 15127 

7.14 + 0.83 

1.43 ± 0.52 

NCC 

RXJ2011.3-5725 

0.279 

24.0 

4995 

4.05 ± 0.49 

0.73 ± 0.27 

CC 

A1758 

0.279 

154.7 

7710 13997 15538 15540 

8.64 + 0.57 

1.36 ± 0.50 

NCC 

A697 

0.282 

19.5 

4217 

11.20 + 1.49 

0.31 ± 0.11 

NCC 


© 0000 RAS, MNRAS 000, 000-000 


The galaxy cluster outskirts probed by Chandra 23 


Table 1. Continued. 


Cluster 

z 

^exp 

ks 

ID 

Tew 

(keV) 

w 

(xlO -2 ) 

CC/NCC 

SDSSJ0922+0345 

0.284 

9.9 

11768 

9.55 ± 2.32 

2.83 + 1.03 

NCC 

RXJ0232.2-4420 

0.284 

22.6 

4993 

9.02 ± 1.18 

2.14 + 0.78 

CC 

A1703 

0.284 

29.1 

15123 

5.95 + 0.33 

0.43 + 0.16 

NCC 

RXJ0437.1+0043 

0.285 

42.5 

11729 7900 

8.10 + 0.85 

0.61 + 0.22 

CC 

ZWCL3146 

0.290 

82.0 

909 9371 

9.98 ± 0.49 

0.46 + 0.17 

CC 

Z7215 

0.292 

13.0 

7899 

6.47 + 1.18 

2.08 ± 0.76 

NCC 

A2813 

0.292 

19.9 

9409 

6.45 + 0.67 

0.74 + 0.27 

NCC 

A2537 

0.295 

38.3 

9372 

6.81 + 0.72 

0.81 + 0.29 

CC 

ACOS0520 

0.295 

31.0 

9331 15099 

8.79 + 1.41 

1.86 + 0.68 

NCC 

1E0657 

0.296 

562.4 

3184 5356 5361 554 4984 4985 4986 5355 

5357 5358 

12.73 + 0.26 

3.10 + 1.13 

NCC 

ACT-CLJ0707-5522 

0.296 

19.8 

12271 

4.77 ± 0.44 

2.52 + 0.92 

NCC 

A781 

0.298 

44.8 

534 15128 

7.24 + 0.71 

3.67+1.34 

NCC 

G292.51+21.98 

0.300 

42.7 

15134 

7.98 + 0.81 

0.83 + 0.30 

NCC 

AS0295 

0.301 

19.8 

12260 

7.72 ± 0.86 

2.08 ± 0.76 

NCC 

MS1008.1-1224 

0.301 

51.2 

926 7711 

7.49 + 0.53 

3.74 + 1.37 

NCC 

A2552 

0.302 

36.3 

11730 3288 

9.30 + 1.07 

0.53 + 0.19 

NCC 

MACS J2245.0+2637 0.304 

14.2 

3287 

6.93 ± 1.29 

0.64 + 0.23 

CC 

Z5699 

0.306 

29.3 

14015 

6.61 + 0.84 

1.53 + 0.56 

NCC 

MACSJ1131.8-1955 

0.307 

23.9 

3276 15300 

10.51 ± 1.13 

3.60 + 1.31 

NCC 

A2744 

0.308 

100.4 

7712 7915 8477 8557 

9.33 + 0.54 

1.82 + 0.66 

NCC 

MACS J0242.5-2132 

0.314 

11.6 

3266 

5.87 + 1.11 

0.38 + 0.14 

CC 

RXCJ2003.5-2323 

0.317 

49.2 

7916 

8.05 + 0.59 

5.61 ± 2.05 

NCC 

MACS J1427.6-2521 

0.318 

28.4 

9373 

6.84 ± 0.89 

0.36 + 0.13 

CC 

A1995 

0.319 

153.3 

7021 7713 7022 7023 

6.79 + 0.55 

0.92 ± 0.34 

NCC 

SPT-CLJ2355-5055 

0.320 

11.5 

11746 

4.90 ± 0.65 

1.05 + 0.38 

CC 

A1351 

0.322 

32.7 

15136 

9.64 ± 1.09 

3.49 + 1.28 

NCC 

MACS J0257.6-2209 

0.322 

20.5 

3267 

6.86 + 1.13 

0.65 + 0.24 

CC 

MACS J0308.9+2645 

0.324 

22.7 

3268 

9.59 + 1.11 

0.86 + 0.31 

CC 

MACSJ2229.7-2755 

0.324 

30.2 

3286 9374 

5.14 + 0.57 

0.18 + 0.06 

CC 

MACS J2135.2-0102 

0.325 

26.7 

11710 

9.31 ± 1.77 

1.39 + 0.51 

NCC 

MACS J2049.9-3217 

0.325 

23.5 

3283 

7.62 ± 1.07 

1.70 + 0.62 

NCC 

ZWCL1358+6245 

0.328 

6.7 

7714 

6.21 ± 1.24 

0.55 + 0.20 

CC 

MACS J0712.3+5931 

0.328 

25.7 

11709 

6.39 + 1.02 

1.32 ± 0.48 

CC 

SPT-CLJ0236-4938 

0.334 

39.6 

12266 

5.40 + 0.80 

1.76 + 0.64 

NCC 

MACSJ0520.7-1328 

0.340 

18.5 

3272 

8.07 + 1.23 

0.51 + 0.19 

CC 

SPT-CLJ2031-4037 

0.342 

9.9 

13517 

6.51 + 0.85 

1.24 + 0.45 

CC 

ACT-CLJ0217-5245 

0.343 

19.8 

12269 

6.37 + 0.91 

3.86 + 1.41 

NCC 

RXJ1532.9+3021 

0.345 

10.0 

1665 

6.47 + 0.99 

0.14 + 0.05 

CC 

SPT-CLJ0040-4407 

0.350 

8.0 

13395 

7.61 + 0.93 

0.87 + 0.32 

NCC 

MACS J1931.8-2634 

0.352 

112.5 

3282 9382 

7.92 ± 0.45 

0.17 + 0.06 

CC 

MACS Jll 15.8+0129 

0.352 

52.9 

3275 9375 

9.14 + 1.06 

0.36 + 0.13 

CC 

A1063S 

0.354 

26.7 

4966 

11.54 ± 1.11 

1.43 + 0.52 

CC 

RBS0797 

0.354 

11.7 

2202 

9.94 + 1.29 

0.21 + 0.08 

CC 

MACSJ0404.6+1109 

0.355 

21.3 

3269 

6.59 + 1.16 

2.07 + 0.75 

NCC 

SPT-CLJ2325-4111 

0.358 

8.9 

13405 

6.50 + 1.17 

3.23 + 1.18 

NCC 

SPT-CLJ0348-4515 

0.358 

12.6 

13465 

5.48 + 0.96 

1.57 + 0.57 

NCC 

MACSJ0011.7-1523 

0.360 

58.9 

3261 6105 

6.21 ± 0.49 

0.53 + 0.19 

CC 

MACS J0150.3-1005 

0.363 

26.8 

11711 

6.75 + 0.85 

0.39 + 0.14 

CC 

MACS J0035.4-2015 

0.364 

21.4 

3262 

8.43 + 0.56 

0.51 + 0.18 

CC 

RXJ0027.6+2616 

0.367 

31.4 

3249 14012 

7.24 ± 1.34 

1.05 + 0.38 

NCC 

CLJ0318-0302 

0.370 

14.3 

5775 

4.25 ± 0.45 

3.44 + 1.26 

NCC 

SPT-CLJ0411-4819 

0.370 

8.0 

13396 

5.63 + 1.01 

0.99 + 0.36 

CC 

ZWCL1953 

0.380 

28.1 

1659 7716 

8.35 ± 1.11 

1.54 + 0.56 

NCC 

MACS J0949.8+1708 

0.384 

14.3 

3274 

9.52 ± 1.87 

1.50 + 0.55 

NCC 

MACS J1731.6+2252 

0.389 

20.5 

3281 

8.97 + 1.74 

1.87 + 0.68 

NCC 

RXJ1720.2+3536 

0.391 

63.2 

3280 6107 7225 7718 

6.50 ± 0.43 

0.60 + 0.22 

CC 

SPT-CLJ0304-4921 

0.392 

20.8 

12265 

6.93 ± 1.25 

0.68 + 0.25 

CC 

MACS J1354.6+7715 

0.396 

32.6 

11754 

6.86 + 1.11 

1.15 + 0.42 

CC 

MACSJ0416.1-2403 

0.397 

15.8 

10446 

8.14 + 1.23 

0.59 + 0.22 

NCC 

MACSJ0429.6-0253 0.399 
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22.9 

3271 

8.55 ± 1.49 

0.47 + 0.17 

CC 





24 Morandi et al. 

Table 1. Continued. 


Cluster 

z 

^exp 

ks 

ID 

Pew 

(keV) 

w 

(xlO -2 ) 

CC/NCC 

MACS J0159.8-0849 

0.406 

72.5 

3265 6106 9376 

8.85 ± 0.87 

0.66 

±0.24 

CC 

SPT-CLJ0013-4906 

0.406 

14.1 

13462 

6.12 ± 1.09 

1.14 

±0.42 

NCC 

CLJ1213+0253 

0.409 

18.4 

4934 

4.29 + 0.88 

1.53 

±0.56 

NCC 

MACSJ2228+2036 

0.412 

19.6 

3285 

9.20 + 1.45 

2.63 

±0.96 

NCC 

SPT-CLJ0234-5831 

0.415 

10.0 

13403 

8.11 ± 1.05 

0.49 

±0.18 

CC 

SLJ1602.8+4335 

0.415 

42.3 

12308 

5.58 + 1.61 

0.66 

±0.24 

NCC 

SPT-CLJ0555-6405 

0.420 

10.9 

13404 

8.65 + 2.04 

2.67 

±0.97 

NCC 

MACS J2046.0-3430 

0.420 

39.2 

9377 

5.73 + 0.71 

0.24 

±0.09 

CC 

SPT-CLJ0252-4824 

0.421 

30.8 

13494 

6.74 + 0.90 

5.59 

±2.04 

NCC 

SPT-CLJ0438-5419 

0.422 

19.8 

12259 

13.17 + 3.18 

1.51 

±0.55 

CC 

SPT-CLJ0551-5709 

0.423 

36.1 

11743 11871 

5.03 + 0.88 

2.47 

±0.90 

NCC 

MS0302.7+1658 

0.424 

10.0 

525 

3.94 + 0.50 

2.07 

±0.75 

CC 

MS1621.5+2640 

0.426 

29.6 

546 

7.18 + 1.38 

1.70 

±0.62 

NCC 

SPT-CLJ2135-5726 

0.427 

17.1 

13463 

6.93 ± 1.18 

0.77 

±0.28 

NCC 

CLJ1216+2633 

0.428 

18.5 

4931 

5.90 + 1.80 

1.15 

±0.42 

NCC 

MACS J0358.8-2955 

0.428 

57.8 

11719 12300 13194 

8.00 + 0.72 

1.88 

±0.69 

CC 

MACS J0451.9+0006 0.430 

10.0 

5815 

7.66 ± 1.04 

2.01 

±0.73 

NCC 

MACS J0553.4-3342 

0.431 

82.9 

12244 5813 

12.68 + 1.26 

2.76 

± 1.01 

NCC 

MACSJ1226.8+2153 0.436 

129.8 

12878 

6.09 + 0.52 

0.38 

±0.14 

NCC 

MACSJ1115.2+5320 0.439 

33.6 

3253 5008 5350 

8.59 + 1.10 

2.24 

±0.82 

NCC 

MACSJ0417.5-1154 

0.440 

89.1 

3270 11759 12010 

10.65 ± 0.54 

2.77 

± 1.01 

CC 

MACS J1206.2-0847 

0.440 

23.2 

3277 

11.41 ± 1.64 

1.31 

±0.48 

CC 

MACS J0913.7+4056 0.442 

76.1 

10445 

6.86 + 0.88 

0.22 

±0.08 

CC 

SPT-CLJ0330-5228 

0.442 

19.4 

893 

4.75 + 0.49 

1.92 

±0.70 

NCC 

MACS J2243.3-0935 

0.447 

20.0 

3260 

9.74 + 1.26 

0.96 

±0.35 

NCC 

MACS J0329.6-0211 

0.450 

69.4 

3257 3582 6108 

5.85 + 0.45 

0.83 

±0.30 

CC 

MACS J1359.1-1929 

0.450 

49.2 

9378 

7.38 + 1.35 

0.76 

±0.28 

CC 

RXJ1347.5-1145 

0.451 

213.8 

3592 13516 13999 14407 

13.63 + 0.69 

0.84 

±0.31 

CC 

MACS J0140.0-0555 

0.451 

29.5 

5013 12243 

8.60 + 1.56 

2.41 

±0.88 

NCC 

RXJ1701+6414 

0.453 

49.3 

547 

6.39 + 0.72 

1.59 

±0.58 

CC 

SPT-CLJ0509-5342 

0.461 

28.8 

9432 

6.49 + 1.21 

0.59 

±0.22 

CC 

MACSJ1621.3+38100.461 

155.0 

3254 3594 6109 6172 9379 10785 

7.41 + 0.67 

0.37 

±0.14 

CC 

CLJ1641+4001 

0.464 

44.5 

3575 

5.03 + 0.68 

4.17 

± 1.52 

NCC 

3C295 

0.464 

87.7 

2254 

5.65 + 0.68 

1.64 

±0.60 

CC 

SPT-CLJ2035-5251 

0.470 

18.4 

13466 

4.90 + 0.84 

0.59 

±0.22 

NCC 

SPT-CLJ0655-5234 

0.470 

20.2 

13486 

8.58 + 1.91 

1.33 

±0.49 

NCC 

CLJ0522-3625 

0.472 

45.3 

4926 5837 

4.94 ± 0.66 

1.42 

±0.52 

NCC 

CLJ0853+5759 

0.475 

42.4 

5765 4925 

7.51 ± 1.85 

1.45 

±0.53 

NCC 

SPT-CLJ2145-5644 

0.480 

14.9 

13398 

7.59 + 1.07 

0.69 

±0.25 

NCC 

ACT-CLJ0215-5212 

0.480 

20.8 

12268 

5.02 ± 1.05 

1.52 

±0.55 

NCC 

SPT-CLJ2233-5339 

0.480 

17.1 

13504 

8.50 + 1.27 

0.89 

±0.32 

NCC 

SPT-CLJ0334-4659 

0.485 

25.7 

13470 

6.44 ± 0.63 

1.13 

±0.41 

CC 

MACS J1824.3+4309 0.487 

14.9 

3255 

4.74 ± 0.97 

3.99 

± 1.46 

CC 

CLJ0926+1242 

0.489 

50.0 

4929 5838 

4.79 ± 0.56 

0.56 

±0.20 

NCC 

MACSJ1427.3+4408 0.490 

40.8 

9380 9808 

8.94 + 1.05 

0.38 

±0.14 

CC 

MACSJ1311.0-0310 

0.494 

114.6 

7721 3258 6110 9381 

6.36 ± 0.58 

0.48 

±0.18 

CC 

SPT-CLJ0200-4852 

0.498 

23.5 

13487 

7.77 + 1.36 

1.82 

±0.67 

NCC 

CLJ0030+2618 

0.500 

16.1 

5762 

5.03 + 0.97 

1.22 

±0.44 

NCC 

MACS J2214.9-1359 

0.503 

32.9 

3259 5011 

10.04 + 1.37 

0.91 

±0.33 

NCC 

MACSJ0257.1-2325 

0.505 

38.1 

1654 3581 

10.41 ± 1.62 

0.70 

±0.26 

CC 

RXJ1525+0958 

0.516 

49.1 

1664 

5.18 + 0.44 

1.52 

±0.55 

NCC 

SPT-CLJ0304-4401 

0.520 

14.9 

13402 

7.43 ± 1.42 

0.95 

±0.35 

NCC 

SPT-CLJ0346-5439 

0.530 

17.8 

12270 

5.15 + 0.69 

1.23 

±0.45 

NCC 

ACT-CLJ0346-5438 

0.530 

16.1 

13155 

4.85 + 0.67 

0.99 

±0.36 

NCC 

X0916+2950 

0.530 

38.3 

12913 

2.81 + 0.76 

1.82 

±0.67 

NCC 

SPT-CLJ2306-6505 

0.530 

25.3 

13503 

5.62 + 0.95 

0.93 

±0.34 

NCC 

MS0016.9+1609 

0.541 

67.4 

520 

9.26 ± 0.77 

0.46 

±0.17 

NCC 

MACS J1149.5+2223 0.544 

38.3 

1656 3589 

8.76 ± 1.05 

1.48 

±0.54 

NCC 

RXJ 1423.8+2404 

0.545 

18.5 

1657 

6.81 + 0.81 

0.37 

±0.13 

CC 

CLJ1354-0221 

0.546 

54.3 

5835 4932 

5.78 + 0.98 

1.64 

±0.60 

NCC 
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Table 1. Continued. 


Cluster 

z 

^exp 

ks 

ID 

Pew 

(keV) 

w 

(xl0“ 2 ) 

CC/NCC 

SPT-CLJ2335-4544 

0.547 

20.7 

13496 

7.67+1.29 

1.27 + 0.46 

NCC 

MACS J0717.5+3745 0.548 

149.1 

1655 4200 16235 

10.55 + 0.55 

1.60 + 0.59 

NCC 

CLJ1117+1744 

0.548 

63.4 

4933 5836 

4.47 + 0.88 

1.59 ± 0.58 

NCC 

SPT-CLJ0307-5042 

0.550 

39.4 

13476 

7.65 + 1.19 

0.48 + 0.17 

NCC 

SPT-CLJ0232-5257 

0.556 

19.8 

12263 

6.87+1.19 

3.54 + 1.29 

NCC 

RXJ1121+2327 

0.562 

69.2 

1660 

4.24 + 0.33 

1.30 + 0.48 

NCC 

SPT-CLJ0456-5116 

0.562 

49.4 

13474 

5.98 + 0.95 

1.15 ± 0.42 

NCC 

SPT-CLJ2148-6116 

0.571 

36.0 

13488 

7.34 + 1.19 

1.75 ± 0.64 

NCC 

CL J0848.7+4456 

0.574 

184.8 

1708 927 

2.72 ± 0.35 

1.67 + 0.61 

NCC 

SPT-CLJ2331-5051 

0.576 

28.7 

9333 

4.80 + 0.86 

0.54 ± 0.20 

CC 

CLJ0216-1747 

0.578 

63.1 

6393 5760 

4.24 + 0.71 

1.52 ± 0.55 

NCC 

SPT-CLJ0307-6225 

0.579 

24.7 

12191 

5.11 + 0.79 

2.21 ± 0.81 

NCC 

SPT-CLJ2245-6206 

0.580 

29.3 

13499 

5.95 + 1.02 

5.72 ± 2.56 

NCC 

MS2053.7-0449 

0.583 

88.6 

1667 551 

4.57 + 0.68 

1.01 ± 0.37 

NCC 

RXJ0647.7+7015 

0.584 

39.3 

3196 3584 

12.57 + 2.10 

1.10 + 0.40 

NCC 

MACS J0025.4-1222 

0.585 

156.9 

10413 10786 10797 3251 5010 

7.94 ± 0.54 

1.34 + 0.49 

NCC 

CLJ0956+4107 

0.587 

58.4 

5759 4930 5294 

5.23 + 0.85 

2.13 ± 0.78 

NCC 

MACS J2129.4-0741 

0.594 

37.8 

3595 3199 

8.11 + 0.89 

2.84 + 1.04 

NCC 

SPT-CLJ2232-5959 

0.594 

31.9 

13502 

5.20 + 0.55 

0.65 ± 0.24 

CC 

SPT-CLJ2344-4243 

0.595 

11.9 

13401 

13.68 + 4.24 

0.21 ± 0.08 

CC 

SPT-CLJ0033-6326 

0.597 

21.1 

13483 

7.60 ± 2.30 

1.65 ± 0.60 

NCC 

CLJ1120+4318 

0.600 

19.8 

5771 

4.80 + 0.52 

1.14 + 0.42 

NCC 

SPT-CLJ0559-5249 

0.609 

108.4 

12264 13116 13117 

5.77 + 0.76 

1.84 + 0.67 

NCC 

CLJ1334+5031 

0.620 

19.5 

5772 

4.99 + 1.35 

0.77 + 0.28 

NCC 

SPT-CLJ0417-4748 

0.620 

21.8 

13397 

6.55 + 0.78 

0.64 ± 0.23 

CC 

SPT-CLJ0123-4821 

0.620 

71.1 

13491 

5.78 + 0.89 

1.52 ± 0.56 

NCC 

SPT-CLJ0256-5617 

0.630 

47.4 

13481 14448 

7.77 + 0.99 

2.44 ± 0.89 

NCC 

SPT-CLJ0426-5455 

0.630 

32.4 

13472 

5.50 + 1.04 

2.02 ± 0.74 

NCC 

CLJ0542.8-4100 

0.634 

50.4 

914 

7.37+1.37 

2.40 + 0.88 

NCC 

SPT-CLJ0243-5930 

0.650 

46.9 

13484 15573 

7.69 ± 0.88 

1.41 ± 0.51 

NCC 

SPT-CLJ2218-4519 

0.650 

34.8 

13501 

8.45 + 1.53 

2.07 + 0.76 

NCC 

SPT-CLJ2222-4834 

0.652 

32.5 

13497 

6.02 + 0.87 

1.16 ± 0.43 

CC 

SPT-CLJ0212-4657 

0.655 

28.0 

13464 

3.84 + 0.61 

3.50 + 1.43 

NCC 

SPT-CLJ0352-5647 

0.660 

44.4 

13490 15571 

5.76 + 0.87 

0.94 ± 0.34 

NCC 

SPT-CLJ0616-5227 

0.684 

38.6 

12261 13127 

6.81 + 0.95 

1.96 ± 0.72 

NCC 

MACS J0744.8+3927 0.686 

87.3 

3197 3585 6111 

8.90 + 0.97 

1.36 ± 0.50 

CC 

RXJ1757 

0.690 

46.1 

10443 11999 

4.12 + 1.14 

1.15 ± 0.42 

NCC 

C11324+3059 

0.696 

48.4 

9403 9840 

4.13 + 2.16 

2.21 ± 0.81 

NCC 

RCS2327.4-0204 

0.700 

141.0 

14025 14361 

10.60 + 0.92 

0.29 ± 0.11 

CC 

RXJ1221+4918 

0.700 

78.3 

1662 

8.10 + 1.25 

2.18 ± 0.79 

NCC 

SPT-CLJ0000-5748 

0.702 

30.1 

9335 

6.54 + 0.76 

0.92 ± 0.34 

CC 

SPT-CLJ0310-4647 

0.709 

36.7 

13492 

7.21 ± 1.38 

1.05 ± 0.38 

NCC 

SPT-CLJ0102-4603 

0.720 

60.7 

13485 

4.56 + 0.80 

4.56 ± 1.66 

NCC 

CL J2302.8+0844 

0.722 

104.7 

918 

6.89 + 1.11 

5.63 ± 2.06 

NCC 

SPT-CLJ2043-5035 

0.723 

79.4 

13478 

6.72 ± 0.74 

0.45 ± 0.17 

CC 

CLJ1113.1-2615 

0.725 

102.3 

915 

5.22 + 0.66 

2.02 ± 0.74 

NCC 

SPT-CLJ0142-5032 

0.730 

29.1 

13467 

8.92 + 2.22 

0.89 ± 0.32 

NCC 

SPT-CLJ2301-4023 

0.730 

57.7 

13505 

5.95 + 0.90 

0.89 ± 0.32 

NCC 

SPT-CLJ2352-4657 

0.730 

79.4 

13506 

8.03 + 1.73 

3.77 + 1.38 

NCC 

SPT-CLJ0406-4805 

0.737 

26.0 

13477 

6.50 + 1.27 

3.83 ± 1.40 

NCC 

SPT-CLJ0324-6236 

0.740 

54.2 

12181 13137 13213 

7.29 + 1.20 

1.42 ± 0.52 

NCC 

SPT-CLJ2259-6057 

0.750 

64.0 

13498 

6.32 + 0.87 

1.96 ± 0.72 

NCC 

SPT-CLJ0014-4952 

0.752 

55.0 

13471 

7.17 + 0.85 

6.79 ± 2.48 

NCC 

C11324+3011 

0.755 

50.2 

9404 9836 

4.69 + 2.30 

1.61 ± 0.59 

NCC 

SPT-CLJ0528-5300 

0.768 

124.0 

11747 11874 12092 13126 9341 10862 11996 

5.26 + 0.78 

0.78 ± 0.28 

NCC 

SPT-CLJ2337-5942 

0.775 

19.7 

11859 

9.12 + 1.52 

0.77 + 0.28 

NCC 

SPT-CLJ2359-5009 

0.775 

131.1 

9334 11742 11864 11997 

6.45 ± 1.08 

1.06 ± 0.39 

NCC 

RCS2318+0034 

0.780 

149.5 

11718 

7.89 + 1.74 

1.06 ± 0.39 

NCC 

SPT-CLJ0449-4901 

0.790 

51.0 

13473 

8.21 ± 1.43 

2.91 ± 1.06 

NCC 

SPT-CLJ0441-4855 0.790 69.2 
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13475 14371 14372 

7.09 ± 0.92 

0.92 ± 0.34 

NCC 
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Table 1 . Continued. 


Cluster 2 

^exp 

ks 

ID 

Tew 

(keV) 

w 

(xlO -2 ) 

CC/NCC 

RXJ1350.0+6007 0.804 

57.8 

2229 

4.79 ± 0.79 

2.63 ± 0.96 

NCC 

RXJ1317.4+2911 0.805 

110.3 

2228 

4.72 + 1.56 

3.49 ± 1.27 

NCC 

RXJ1716.9+6708 0.813 

51.2 

548 

6.93 + 1.08 

1.18 ± 0.43 

NCC 

RXJ1821.6+6827 0.820 

49.6 

10444 10924 

4.70 + 0.81 

1.85 ± 0.68 

NCC 

SPT-CLJ0058-6145 0.830 

50.6 

13479 

5.76 ± 1.40 

1.06 ± 0.39 

NCC 

SPT-CL J2258-4044 0.830 

54.4 

13495 

7.30 + 1.07 

1.80 + 0.66 

NCC 

CLJ0152.7-1357 0.831 

35.4 

913 

4.85 + 0.95 

2.70 + 0.99 

NCC 

SPT-CL J0102-4915 0.870 

59.2 

12258 

13.17 + 1.35 

7.10 + 2.59 

CC 

SPT-CL J0533-5005 0.881 

71.7 

11748 12001 12002 

6.11 + 2.05 

1.72 ± 0.63 

NCC 

CL J 1226.9+3332 0.890 

63.4 

3180 5014 

11.68 + 1.73 

1.09 ± 0.40 

CC 

01604 0.898 

94.5 

6932 6933 7343 

4.30 + 2.07 

1.83 + 0.67 

NCC 

SPT-CL J2034-5936 0.920 

58.7 

12182 

7.92 + 1.35 

1.13 ± 0.41 

NCC 

SPT-CL J2146-4632 0.932 

81.0 

13469 

5.24 + 0.62 

1.58 + 0.58 

NCC 

PLCKG266.6-27.3 0.940 

240.5 

14017 14018 14349 14350 14351 14437 15572 15574 

15579 15582 15588 15589 

10.32 + 0.50 

0.77 + 0.28 

CC 

SPT-CL J2345-6405 0.940 

65.2 

13500 

6.59 + 0.92 

1.66 ± 0.61 

NCC 

SPT-CLJ2341-51191.000 

80.0 

11799 9345 

6.32 + 0.85 

1.47 + 0.54 

NCC 

SPT-CL J0037-50471.026 

72.5 

13493 

3.22 + 1.05 

3.86 ± 1.41 

NCC 

CLJ1415.1+3612 1.030 

88.6 

4163 

6.74 + 0.69 

1.05 ± 0.38 

NCC 

SPT-CLJ0151-59541.030 

49.4 

13480 

4.02 + 0.62 

1.53 ± 0.56 

NCC 

SPT-CL J0546-53451.066 

69.5 

9332 9336 10851 10864 11739 

8.54 + 1.11 

1.02 ± 0.37 

NCC 

SPT-CL J2342-54111.075 

174.0 

11741 11870 12014 12091 

4.70 + 0.61 

1.43 + 0.52 

CC 

RXJ0910+5422 1.110 

169.8 

2452 2227 

4.58 + 1.04 

1.54 + 0.56 

NCC 

SPT-CL J2106-58441.132 

24.7 

12180 

8.80 + 1.44 

1.45 ± 0.53 

NCC 

SPT-CL J0446-58491.160 

53.4 

13482 15560 

6.13 + 2.03 

2.44 ± 0.89 

NCC 

SPT-CL J2236-45551.160 

83.0 

13507 15266 

7.44 ± 1.39 

2.85 ± 1.04 

NCC 

SPT-CL J0156-55411.220 

77.7 

13489 

7.81 ± 1.23 

0.75 ± 0.28 

NCC 

RDCSJ1252-29 1.237 

187.3 

4198 4403 

6.50 + 1.08 

1.52 ± 0.56 

NCC 
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